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We study the radiative charmless B^ — > V{A) r ) decays in perturbative QCD (pQCD) 
approach to the leading order in a s (here V and A denotes vector mesons and two kinds of 
axial- vector mesons: 3 Pi and 1 P\ states, respectively). Our predictions of branching ratios 
arc consistent with the current available experimental data. We update all — > V form 
factors and give the predictions for B — > A form factors using the recent hadronic inputs. In 
addition to the dominant factorizable spectator diagrams, which is form factor like, we also 
calculate the so-called "power suppressed" annihilation type diagrams, the gluonic penguin, 
charming penguin, and two photon diagrams. These diagrams give the main contributions 
to direct CP asymmetries, mixing-induced CP asymmetry variables, the isospin asymmetry 
and U-spin asymmetry variables. Unlike the branching ratios, these ratios or obscrvables 
possess higher theoretical precision in our pQCD calculation, since they do not depend on 
the input hadronic parameters too much. Most of the results still need experimental tests 
in the on-going and forthcoming experiments. 
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I. INTRODUCTION 

Present studies on B decays are mainly concentrated on the precise test of standard model (SM) 
and the search for signals of possible new physics. Radiative processes b — > 57 and b — > dry are 
among the most ideal probes and thus have received considerable efforts In the standard model, 
these processes are induced by the flavor-changing-neutral-current transitions which are purely loop 
effects. Such decays are rare and the measurement of parameters in these channels, especially the 
CP asymmetry variables, may shed light on detailed information on the flavor structure of the 
electroweak interactions. This predictive power relies on the accuracy of both the experimental 
side and the theoretical side. Thanks to the technique of operator-product-expansion, remarkable 
theoretical progress has been made in the SM to the next-to-next-to-leading order accuracy [2|. 
Compared with experimental results [jj, the theoretical prediction is consistent with but 1.2a below 
the experimental data 4|. This consistence could certainly provide a rather stringent constraint 
on non-standard model scenarios. 

Compared with inclusive processes, the exclusive processes B — ► V7 is more tractable on the 
experimental side, but more difficult on the theoretical side. Theoretical predictions are often 
hampered by our ability to calculate the decay amplitude (V ~f\Oi\B) , where Oi is a magnetic 
moment or a four-quark operator. We have to use some non-perturbative hadronic quantities to 
describe the bound state effects. In the heavy-quark limit, the non-perturbative contributions can 
be organized in a universal and channel-independent manner. Factorization analysis, the separation 
of the short-distance and long-distance dynamics, can give many important predictions. 

The dominant contribution to the radiative decay amplitudes is proportional to the transition 
form factor. Different treatments on the dynamics in form factor F B ^ V result in different explicit 

n 

approaches. There are many approaches such as Lattice QCD (LQCD) [5j, light-front quark model 
(LFQM) and light-cone sum rules (LCSR) [§]. Recently three commonly- accepted approaches 



are developed to study exclusive B decays: QCD factorization (QCDF) [9(, soft collinear effective 
theory (SCET) [u| and perturbative QCD approach 11 j. In pQCD approach, the recoiling meson 



moves on the light-cone and a large momentum transfer is required. Keeping quarks' transverse 
momentum, pQCD approach is free of endpoint divergence and the Sudakov formalism makes it 
more self-consistent. It has been successfully applied to various decay channels and quite recently 
this approach is accessing to next-to-leading order accuracy jl^ . A bigger advantage is that we can 
really do the form factor calculation and the quantitative annihilation type diagram calculation in 
this approach. The importance of annihilation diagrams are already tested in the predictions of 



direct CP asymmetries in B° — > 7r + 7r~, i^ + 7r~ decays [ill. [13I] and in the explanation of £? — > 
polarization problem Q, [2]. These "power suppressed" contributions are the main source of 
isospin symmetry breaking (or SU(3) breaking) effects in radiative decays B — > V7 [l6. fiyl 
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Some of charmless B — > V7 decay channels have been studied in pQCD approach 
separately. According to the transition at the quark level, all the 10 decay channels can be divided 
into three different groups: b — > s, b —> d and purely annihilation type. The two B — > i^*7 modes 
(6 — > s process) have been well measured experimentally. The agreement of pQCD approach results 
and experimental data is very encouraging. In this paper, we study all those channels including 
the corresponding B s decays in a comprehensive way. The input hadronic parameters will be 
chosen the same as that in Ref. [ijj]. We also updated all the B — > V decay form factors in pQCD 



approach [20J. Despite branching ratios and CP asymmetry parameters which heavily depend on 
the input parameters, we also study some ratios characterizing the isospin and SU(3) breaking 
effects, where most of the uncertainties cancel. Experimental measurements of these quantities 
will prove a good test of our theory, since different method gives different prediction. With the 
ongoing B factories BaBar and Belle, the £>-Physics program on CDF and the onset of the LHC 
experiments, as well as the Super B-factories being contemplated for the future, we expect a wealth 
of data involving these decays. 

Although experimentalists have already measured one channel of the B — > A"f decays [211 ] . there 



are not many discussions on the theoretical side. The pQCD study on B — > V7 can be straight- 
forward extended to radiative processes involving higher resonants such as ETi(1270), i^i(1400). 
In the quark model, the quantum numbers J PC for the orbitally excited axial-vector mesons are 
1 ++ or 1 H , depending on different spins of the two quarks. In SU(3) limit, these mesons can not 
mix with each other; but since the s quark is heavier than u, d quarks, Xi(1270) and Xi(1400) 
are not purely l 3 Pi or l x Pi states. These two mesons are believed to be mixtures of K\a and 
Kib, where K\a and K\b are 3 Pi and l P\ states, respectively. In general, the mixing angle can 
be determined by experimental data. But unfortunately, there is not too much data on the mixing 
of these mesons which leaves the mixing angle much free. Analogous to 77 and r]' , the flavor-singlet 
and flavor-octet axial-vector meson can also mix with each other. Using those hadronic parameters 
determined in B — > V"f decays, the production of A'y in B decays could provide a unique insight 
to these mysterious axial- vector mesons. 

This paper is organized as follows: In section II, we briefly review pQCD approach with the 
operator basis used subsequently. Some input quantities which enter pQCD approach, wave func- 
tion of the £>-meson, distribution amplitudes for light vector mesons, and for light axial-vector 
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mesons and input values of the various mesonic decay constants, are also given here. In section HTT1 
we give the factorization formulae and the numerical results for B — »• V and B — > A form fac- 
tors. Section IIVI contains the calculation of B — ► Vj decays, making explicit the contributions 
from the electromagnetic diploe operator, the chromo-magnetic moment operator, some higher or- 
der {0{a s )) corrections from tree operators, the contribution from two-photon diagrams, the tree 
operator annihilation diagrams and penguin operator annihilation diagrams. Numerical results 
for the charge-conjugated averages of decay branching ratios are given in comparison with the 
corresponding numerical results obtained in QCDF approach and SCET, as well as the available 
experimental data. We also give direct CP-asymmetries, time-dependent CP asymmetries Sf and 
observables Hf (for B s system) in the time-dependent decay rates in this section. In section [Vj we 
study B — > A"f decays by predicting branching ratios and direct CP asymmetries. Our summary is 
given in the last section. Some functions are relegated to the appendix: appendix A contains the 
various functions that enter the factorization formulae in the pQCD approach; appendix B and C 
give the analytic formulae for the B — > Vj and B — > Aj decays used in numerical calculations, 
respectively. 

II. FORMALISM OF PQCD APPROACH 
A. Notations and conventions 



We specify the weak effective Hamiltonian which describes b — > D {D = d,s) transitions 221 ] : 
H *ff = 7§{ E V qb V* D [C 1 { l i)0\{ l x) + 02^)01^) 

V V -y n, r . 



q=u,c 

10,77 



iu,i7,o(/ 

-W[ £ c i ( JJl )o i (ji)]\+n.c., (i) 

i=3 ^ 

where V q b(D) an d Vtb(D) are Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. Functions Oi 
{i = 1, 10, 77, 8g) are local four-quark operators or the moment type operators: 

• current-current (tree) operators 

0\ = (q a bf 3 ) v ^ A (D l3 q a ) v _ A , 0\ = (q a b a ) v _ A (Di3qf3)v-A, (2) 

• QCD penguin operators 

O3 = (D a ba)v-A^2(cfpq'p)v-A, 4 = (Dpb a ) V -A^2(q'a<l'p)v-A, (3) 

q> q> 

O5 = {D a b a ) v _ A ^{q' p q' p ) v+A , 6 = {D p b a ) V - A ^2(<f a tfp)v+A, (4) 

q> q> 
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electro-weak penguin operators 



2 _ 3 _ 

7 = -(D a b a )v-A^eq>{ctp<jp)v+A, 08 = -(Df 3 b a ) V - A ^2e g ,(q' a q' f3 ) v+ A, 



q' 



q' 



(5) 
(6) 



magnetic moment operators 



77 = -^D a a^ v {m D P L + m b P R )b a F^ v , O s 



Air 2 



D a a^(m D P L + m h Pa)T^gbpG%(7) 



where a and (3 are color indices and q' are the active quarks at the scale m&, i.e. q' = (u, d, s, c, b). 
The left handed current is defined as (^q'^v-A = 9^7^(1 — 7s)<70 an d the right handed current 
(q' a q'p)v+A = ^,71/(1 + 75)9/3- The projection operators are defined as Pl = (J^ — 7s)/2 and 
Pr — (1 +75)/2. The combinations aj of Wilson coefficients are defined as usual 231 ]: 



ai = C 2 + Ci/3, a 2 = Ci + C 2 /3, a 3 = C 3 + C 4 /3, 04 = C 4 + C 3 /3, a 5 = C 5 + C 6 /3, 

a 6 = C 6 + C5/3, a 7 = C 7 + C 8 /3, a 8 = C 8 + C 7 /3, o 9 = C 9 + d /3, a w = C w + C 9 /3. (8) 



For the explicit formulae, we will consider B meson decays and use light-cone coordinates to 
describe the momentum: 



,p0 _|_ p3 pO _ p3 
P={P ,P ,PT) = [ — , —>PT 



* V2 ' V2 '"V- (9) 
where pr = ip l ,P 2 )- I n the P< meson rest frame, momenta of B meson, vector (axial- vector) meson 
and the photon are chosen as: 



Pi = -1(1,1, Or), P 2 = ^|(0,l,0r), P 7 = -1(1,0, Or), 



(10) 



where the vector (axial- vector) is mainly moving on the minus direction ra_ and the photon is 
moving on the plus direction n+ . Longitudinal momenta fractions of the spectator anti-quarks in 
B and final state meson are chosen as x\ = kf jP^ and x 2 = k^ /P 2 ~ ■ 1 Including the transverse 
components, the momenta of these spectator antiquarks are expressed by: 



fei = (— j=xi,0, h T ), k 2 = (0, ~^= x 2, kzr), 



(11) 



1 One should be cautious that in the discussion of light cone distribution amplitudes (LCDAs) for the vectors and 
axial- vectors, x is defined as the momentum fraction of the positive quark which is different with our definition 
here. 



D 6 ' — o? ~, A ' D 



FIG. 1: Feynman diagrams of the electromagnetic penguin operator 

then the b and D{= s, d) quark momenta are Pb = Pi — k\ and pr> = P2 — &2- For convenience, we 
can define the following useful ratio variables: 

m b m D my m A /10 , 
n = , r D = , r v = , r A = , (12) 

mB m-B n^B n^B 

where mo and m^^) are masses for the d(s) quark and the vector (axial- vector) meson, respec- 
tively. In the calculation of decay amplitudes, we will only keep terms of leading order in ry but we 
will consider the corrections together with kinematic corrections in the phase space as in Eq. (|87p . 

According to the Lorentz structure, decay amplitudes from various operators can be generally 
decomposed into scalar and pseudoscalar components as: 

M = {e* ■ e* v )M s + ie llva p^€^n%n p _M p , (13) 

and where we adopt the convention e 0123 = +1. In order to study mixing- induced CP asymmetries, 
it is convenient to separate different chiralities in the amplitudes. If the emitted photon is left- 
handed, the relationship between the scalar M. and pseudoscalar component M. p is required 

as 

M s = -M p = ^M L , (14) 



while the condition 



M S = M P = ^M R , (15) 



is required for the right-handed photon. 



B. A brief review of pQCD approach 



The basic idea of pQCD approach is that it takes into account the intrinsic transverse momentum 
of valence quarks. The decay amplitude, taking the first diagram in Fig. [JJ as an example, can 
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FIG. 2: 0(a s ) corrections to the hard scattering kernel H . 

be expressed as a convolution of wave functions (f) B , 4>v an d hard scattering kernel Tjj by both 
longitudinal and transverse momenta: 
1 r d 2 k 1T d 2 k 2T 



M= lo dXldX2 J (2tt») (2r») 



<f> B (xi, k 1T ,pi, t)T H (xi,x 2 , h T , k 2 T,t)4> v (x 2 , k 2T ,p 2 ,t). (16) 



Usually it is convenient to compute the amplitude in coordinate space. Through Fourier transfor- 
mation, the above equation can be expressed by: 



M = J dx 1 dx 2 I d hd b 2 4> B (xi,bi,pi,t)T H (x 1 ,x 2 ,bi,b 2 ,t)(f>v(x 2 ,b 2 ,p 2 ,t). 



(17) 



This derivation is mainly concentrated on tree level diagrams, but actually we have to take into 
account some loop effects which can give sizable corrections. The 0(a s ) radiative corrections to 
hard scattering process H are depicted in Fig. [2j In general, individual higher order diagrams may 
suffer from two types of infrared divergences: soft and collinear. Soft divergence comes from the 
region of a loop momentum where all it's momentum components vanish: 



F = (l+,r,l T ) = (A,A,A), 



(18) 



where A is the typical scale for hadronization. Collinear divergence originates from the gluon 
momentum region which is parallel to the massless quark momentum, 



Z" = (l+,r,T T ) ~ (m B ,A 2 /m B , A). 



(19) 



In both cases, the loop integration corresponds to f d l/l 4 ~ log A, thus logarithmic divergences 
are generated. It has been shown order by order in perturbation theory that these divergences 
can be separated from the hard kernel and absorbed into meson wave functions using eikonal 
approximation 



24f | . But when soft and collinear momentum overlap, there will be double logarithm 
divergences in the first two diagrams of Fig. [2j These large double logarithm can be resummed 
into the Sudakov factor whose explicit form is given in Appendix lAl 

Furthermore, there are also another type of double logarithm which comes from the loop correc- 
tion for the weak decay vertex correction. The left diagram in Fig. [T]gives an amplitude proportional 
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to l/(x%xi). In the threshold region with X2 — > [(to be precise, X2 ~ 0(Aqcd /™b))], additional 
soft divergences are associated with the internal quark at higher orders. The QCD loop corrections 
to the electro-weak vertex can produce the double logarithm a s In 2 X2 and resummation of this 
type of double logarithms lead to the Sudakov factor St(x2)- Similarly, resummation of a s In 2 x\ 
due to loop corrections in the other diagram leads to the Sudakov factor St(x\). These double 
logarithm can also be factored out from the hard part and grouped into the quark jet function. 
Resummation of the double logarithms results in the threshold factor 25[. This factor decreases 



faster than any other power of x as x — ► 0, which modifies the behavior in the endpoint region to 
make pQCD approach more self-consistent. For simplicity, this factor has been parameterized in a 
form which is independent on channels, twists and flavors 26]. 



Combing all the elements together, we can get the typical factorization formulae in pQCD 
approach: 

M = j dxidx 2 J d 2 bid 2 b2{^) 2 (pB{xi,bi,pi,t) 

xT H (x 1 ,x 2 ,Q,bi,b2,t)(j) V (x2,b2,P2,t)St(x2)eyLp[-SB(t) - S 2 (t)\. (20) 

C. Wave functions of B mesons 



In order to calculate the analytic formulas of decay amplitudes, we will use light cone wave 
functions <&M,a/3 decomposed in terms of the spin structure. In general, &M,a/3 with Dirac indices 
a,f3 can be decomposed into 16 independent components, l a /3, 75 Q /?, 7^0, (7^75)0/3, cr^g- If the 
considered meson is the B meson, a heavy pseudo-scalar meson, the B meson light-cone matrix 



element can be decomposed 



27 



28, 



m by: 



d*ze ik ^{0\bp(0)D a (z)\B(P B )) 



[Pb + m B )j5 



V2 



-<p B {ki 



where N c = 3 is the color factor, n and v are two light-like vectors: n 2 



(21) 



0. From equation 



(|2ip . one can see that there are two Lorentz structures in the B meson distribution amplitudes. 
They obey the following normalization conditions: 



/ 



d 4 h 



4> B {ki 



Sb 



d A h 



>B{Kl, 



0. 



(22) 



2 v / 2lVc' J (2vr) 4 

In general, one should consider these two Lorentz structures in the calculations of B meson 
decays. However, it is found that the contribution of <f> B is numerically small [2^], thus its con- 
tribution can be safely neglected. With this approximation, we only retain the first term in the 
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square bracket from the full Lorentz structure in Eq. ()2ip : 

<S> B = -^{Fb + msHMh)- (23) 

In the following calculation, we will see that the hard part is always independent of one of the kf 
and/or k^ ■ The B meson wave function is then a function of the variables k^ (or kf) and k\ only, 

MKMt) = J^tMK.KM- (24) 

In the 6-space, the B meson's wave function can be expressed by 

$ B (x, b) = -j^= []P B <y5 + m B75 ] cj) B {x, b), (25) 

where b is the conjugate space coordinate of the transverse momentum kj<- 
In this study, we use the following phenomenological wave function: 



4>b(x, b) = Nbx 2 (1 — x) 2 exp 



(26) 



with Nb the normalization factor. In recent years, a lot of studies have been performed for B^ and 
B^ decays in pQCD approach. The parameter loi, = (0.40 ± 0.05) GeV has been fixed using the 
rich experimental data on B® and B^ meson decays. In the SU(3) symmetry limit, this parameter 
should be the same in B s decays. Considering a small SU(3) breaking, s quark momentum fraction 
should be a little larger than that of u or d quark in the lighter B mesons, since s quark is heavier 
than u or d quark. We will use u>b = (0.50 ± 0.05) GeV in this paper for the B s decays 



D. Light-cone distribution amplitudes of light vector mesons 

Decay constants for vector mesons are defined by: 

(0| 917^92 1 V(p, e)) = fvmve^ (0|gio><7 2 |V(p, e)) = ifvi^Pv - e^). (27) 

The longitudinal decay constants of charged vector mesons can be extracted from the data on 
r~ — ► (p~ , K*~)v T jsfj. Neutral vector meson's longitudinal decay constant can be determined by 
its electronic decay width through V° — > e + e~ 2 and results are given in Table HI Transverse decay 
constants are mainly explored by QCD sum rules [31[ that are also collected in Table HI 



2 There is a recent study on extracting vectors' decay constants from experimental data, which has taken into 
account the effects of p°-ui and u>-(f> mixing [3f[. Since we do not consider the mixing for the decay amplitudes in 
our calculation, we will not use those values for self-consistence. 
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TABLE I: Input values of the decay constants for the vector mesons (in McV) 

fp fj fu fx* fx* ftp f<p 

209 ± 2 165 ± 9 195 ± 3 151 ±9 217 ± 5 185 ± 10 231 ± 4 186 ± 9 



We choose the vector meson momentum P with P 2 = m v , which is mainly on the plus direction. 
The polarization vectors e, satisfying P ■ e = 0, include one longitudinal polarization vector el and 
two transverse polarization vectors €t- The vector meson distribution amplitudes up to twist-3 are 
defined by: 

1 ^ 



(V(P,e* L )\q 2p (z)q la (0)\0) 
(V(P,e* T )\q 2 p(z)q la (0)\0) 



V2WcJo 



dxe ixP - z [m v f L M*)+ <L PtfvW+mvM*)] 



a/3 



-±= Cdxe^ [my f T <P v {x)+ f T Pf v {x) 

y/Zl\ c Jo 

+m v ie^ pul ^e?n p v (J <i)y(x)} al3 , 



(28) 



for the longitudinal polarization and transverse polarization, respectively. Here x is the momentum 
fraction associated with the q 2 quark, n is the moving direction of the vector meson and v is the 
opposite direction. These distribution amplitudes can be related to the ones used in QCD sum 
rules by: 

fv 



<t>v(x) 

4> V {X) 



-4\\{x), 4>y{x) 



fT 



fv d (s 



f T 



<j> v {x) = g<p (x), (f> v (x) = ^ 9 { f (x) . 



(29) 



The twist-2 distribution amplitudes can be expanded in terms of Gegenbauer polynomials C n 
with the coefficients called Gegenbauer moments a n : 



3/2 



<P\\,±(x) = 6x(l - x) 



n=l 



(30) 



where t = 2x — 1. The Gegenbauer moments a)n are main 
sum rules. Here we quote the recent numerical results [321 . 



y determined by the technique of QCD 



33. 



351 ] as 



a\(K*) = 0.03 ±0.02, af(K*) = 0.04 ± 0.03, 

al(p) = o|M = 0.15 ±0.07, a^(p) = a^{u) = 0.14 ± 0.06, 
a\(K*) = 0.11 ±0.09, c£ (K*) = 0.10 ± 0.08, 



4(0) 



0.18 ±0.08, 



a 2 



0.14 ±0.07, 



(31) 
(32) 
(33) 
(34) 
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TABLE II: Input values of the decay constants (absolute values) for the axial- vector mesons (in MeV). The 
transverse decays constants for 1 Pi are evaluated at \i = 1 GeV. 

/in (1260) //i(l 3 Pi) //e(l 3 Pi) /Kia /£(1235) fh 1 {l l P 1 ) fh 8 (l 1 P 1 ) f K XB 

238 ± 10 245 ± 13 239 ± 13 250 ± 13 180 ±8 180 ± 12 190 ± 10 190 ± 10 

where the values are taken at fi = 1 GeV. 

Using equation of motion, two-particle twist-3 distribution amplitudes are related to twist- 2 
LCDAs and three-particle twist-3 distribution amplitudes. But in some B — > VV decays, there 
exists the so-called polarization problem. It has been suggested that using asymptotic LCDAs can 
resolve this problem in pQCD approach. Thus to be self-consistent, we should also use the same 
form to calculate radiative decays. As in Ref. we use the asymptotic form for twist-3 LCDAs: 

hf(x) = 3t 2 , h^\x) = 6x(l-x), (35) 
g { l\x) = 6x(l-x), g { l\x) = \{l + t 2 ). (36) 

E. Light-cone distribution amplitudes of axial-vectors 

Longitudinal and transverse decay constants for axial-vectors are defined by: 

(A(P, e)|g 27M 75<7i|0> = i/im A£ ;, (A(P, e)|g 2 <x^ 7 5gi|0) = /l(e*P„ - €* U P, L ). (37) 

In SU(2) limit, due to G-parity, the longitudinal (transverse) decay constants vanish for the 
non-strange 1 Pi[ 3 Pi] states. This will affect the normalization for the corresponding distribu- 
tion amplitudes which will be discussed in the following. For convenience, we take fzp 1 = f 
[fip (/i = 1 GeV) = /] as the "normalization constant" . The numbers of axial vector meson decay 
constants shown in table [TT] are taken from Ref. [3^] • 

Distribution amplitudes for axial-vectors with quantum numbers J PC = 1 ++ or l" 1 are defined 

by: 

(A(P, e* L )\q 2f3 (z)q la (0)\0) = -J= J dxe**' [m A f L ls<t> A {x)- f L P-%<t>\{x) - m Al ^ A {x)] ap , 

(A(P,e* T )\q 2 p(z)q la (0)\0) = £ dxe^* [m A f T i^\{x)- f T f>^ T A {x) 

-m A ie lwpa ^nPv a <f> a A (x)] aj3 . (38) 

Besides the factor —i-fs from the right hand, axial-vector mesons' distribution amplitudes can be 
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related to the vector ones by making the following replacement: 



<t>v -> 4>a, 4>v -<I>ai <Pv -<Pai 

<$^-<f>l, <f>v^<t>% (39) 

These distribution amplitudes can be related to the ones calculated in QCD sum rules by: 

where we use / as the "normalization" constant for both longitudinal polarized and transversely 
polarized mesons. 

In SU(2) limit, due to G-parity, ^>ii, gj° and g± are symmetric [antisymmetric] under the 
replacement x <-> 1 — x for non-strange l 3 Pi [l 1 ^] states, whereas (f)±, and hft are antisym- 
metric [symmetric]. In the above, we have taken = /ap 1 = / [fip l = fT Pl {^ = 1 GeV) = /], 
thus we have 

(1 3 Pl(P, 6)191(7^9210) = fT Pl %' 3Pl (e*P„ - e*P M ), (41) 
(^^(P.eJIgn^^lO) = i/ 1Pl ajj' 1Pl m 1Pl e*, (42) 

J- 3 Pi II 1 -Pi 

where a ' and a Q are the Gegenbauer zeroth moments. That can give the following normal- 
ization for the distribution amplitudes: 



L 



i 

dx(f>±(x) = ciq 



dxcfrn (x) = ag 



(43) 



for the 1 3 P\ [l^Pi] states. The zeroth Gegenbauer moments a^' Pl and ajj' Pl , characterizing the 
degree of the flavor SU(3) symmetry breaking, are non-zero for only strange mesons. We normalize 
the distribution amplitude 4>\\ of the l 3 Pi states as 

/ dxcj)n(x) = 1. (44) 
Jo 

For convenience, we formally define ajj = 1 for the l 3 Pi states so that we can use Eq. (|43|) as 
the normalization condition. Similarly, we also define Oq- = 1 for l*Pi states so that (ft±(x) has a 
correct normalization. 
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TABLE III: Gegenbauer moments of cf>± and <f>\\ for l 3 Pi and l 1 ^ mesons evaluated in Ref. 37], where the 
values are taken at fx = 1 GeV. 



||,oi(1260) 
"2 




in 


4^ 




-0.02 ±0.02 


-0.04 ±0.03 


-0.07 ±0.04 


-0.05 ±0.03 


0.00 ±0.26 


J_,ai(1260) 
a l 




-L./a 








-1.04 ±0.34 


-1.06 ±0.36 


-1.11 ±0.31 


-1.08 ±0.48 


0.08 ±0.09 


0.02 ±0.20 


||,bi(1235) 
"1 


y-- 










-1.95 ±0.35 


-2.00 ±0.35 


-1.95 ±0.35 


-1.95 ±0.45 


0.14 ±0.15 


0.02 ±0.10 


_L,&i(1235) 


a 2 


a 2 8 




±,KlB 

"1 


0.03 ±0.19 


0.18 ±0.22 


0.14 ±0.22 


-0.02 ±0.22 


0.17 ±0.22 



Up to conformal spin 6, twist- 2 distribution amplitudes for axial- vector mesons can be expanded 



as: 



0||(x) 

4>±{x) 



6xx 
6xx 



ajj + 3a|t + 4 ~ !) 

Oq + 3aft + 4 -(5t 2 - i; 

n 



(45) 
(46) 



wher the Gegenbauer moments are calculated in Refs. [sf], |37( shown in table IIIIl From the 
results in table IIIIl we can see that there are large uncertainties in Gegenbauer moments which 
can inevitably induce large uncertainties to branching ratios and CP asymmetries. We hope the 
uncertainties could be reduced in future studies in order to make more precise predictions. 
As for twist-3 LCDAs, we use the following form: 



+ t 2 ) + - a\ t\ (x) = 6xx(aJ + a\t) 



h^ix) = 3a^t 2 + -a±t(3t 2 - 1), h ( °\x) = 6xx(a^ + aft). 



(47) 
(48) 



14 



III. B -> V AND B -> A FORM FACTORS 



i? — ► V form factors are defined under the conventional form as follows: 

2V(q 2 ) 



(V(P 2 ,e*)\qi fl b\B(P 1 )) 



TUB + my 



(V(P 2 ,e*)\q^ l5 b\B(P 1 )) = 2im v A (q 2 ) e —^-q» + i(m B +m v )A 1 {q 2 ) 



-iA 2 {q 2 )- 



e -g 



(Pi + W 



mjj — m 



me + my 

(V(P 2 ,e*)|^^6|S(P 1 )) = -2iTt(q 2 )e^etP lp P 2a , 
(y{P 2 ,e*)\qa^ imv b\B{P x )) = T 2 (q 2 ) [(m 2 B - m y )e* M - (e* • q)(Pi + P 2 Y 



+W)(e*-«) 



2 2 

m^ — rriy 



(Pi + P2 



(49) 



where q = P\ — P 2 and the relation 2my>lo(0) = (mg + my)j4i(0) — (rne — mv)A 2 (0) is obtained 
in order to cancel the pole at q 2 = 0. 

In pQCD approach, the factorization formulae for these form factors at maximally recoiling 
(q 2 = 0) are expressed by: 



V 



87rCirms(ms + my) / dx\ 
Jo 



dx 2 / bidbib 2 db 2 4>B(xi,bi] 
Jo 

' lT 



y.{ E e (t a )h e (x 1 X2,x 2 ,b l ,b 2 ) 4> v (x 2 ) + (2 + x 2 )r v 4> v (x 2 ) - r v x 2 <f) v (x 2 ) 



+E' e (t' a )h e (xix 2 , x u b 2 , bi)r v }0 v (x 2 ) + 4> v v {x 2 ) 

Aq = &nCFm 2 B I dx\dx 2 I bidb\b 2 db 2 cf)B(xi,bi) 
Jo Jo 

x\ E e (t a )h e (xix 2 ,x 2 ,bi,b 2 ) (1 + x 2 )(f> v (x 2 ) + (1 - 2x 2 )r v ((f>v (^2) + <t>v( x 2)) 



(50) 



+2r v E e (t a )h e (x 1 x 2 , xt,b 2 , bi)4> s v (x 2 ) 

r 1 roo 

Ai = 87rCi?ms(mB — my) / dx\dx 2 / b\db\b 2 db 2 4>B{x\,b\) 

Jo Jo 

x \ E e (t a )h e (xix 2 ,x 2 ,bi,b 2 ) 4> v (x 2 ) + (2 + x 2 )r V (/) v (x 2 ) - r v x 2 cf> v (x 2 ) 



(51) 



A, 



+E' e (t' a )h e (xix 2 ,x 1 ,b 2 ,b 1 )r v <j) v (x 2 ) + 4> v (x 2 ) 



(iris + my)j4i — 2mvAo 



m b — my 



(52) 
(53) 
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8irCFm B / dx\dx 2 / bidbib 2 db 2 4> B (xi,bi) 
Jo Jo 

x \ E e (t a )h e ( Xl X2,x 2 , 61, 6 2 ) [(1 + x 2 )(fv{x2) + (1 - 2x 2 )rv(^v-(a;2) + 4fy{x 2 )) 



+£ , e (t a )/l e (xiX2,Xi,62,6i)ry((/)y(x2) + <^y(x 2 )) 

T3 = T 2 - lQ-nCFm 2 B ry / dx\dx 2 / 61^6162^62^(^1,61) 

io jo 

x<^ E e (t a )h e (xxx 2 ,x 2 ,b 1 ,b 2 ) (f>v(x 2 ) + 2r V (j) t v (x 2 ) + rya^y (^2) - </>v(x 2 )) 



(54) 



+2ry J Eg(i' a )/i e (a;i^2, £1, 6 2 , h)(p s v (x 2 ) >, 



(55) 

with Cp = 4/3. The definitions of functions Ei, the factorization scales ti and hard functions hi, 
are given in Appendix lAl 

With the above expressions for form factors, we obtain the numerical results that are collected 
in table IIVI When evaluating the form factor Ao, we used 2my^4o(0) = (m B + my)j4i(0) — (m B — 
mv)A2(0) and assume Ai and A2 are linearly correlated to estimate the uncertainties. The first 
error comes from decay constants of meson and shape parameters uj^, while the second one is 
from hard scale t and A QC d- In the calculation, f B = (0.19 ± 0.02) GeV, uj b = (0.40 ± 0.05) GeV 
(for B ± and B° d mesons) and f Bs = (0.23±0.02) GeV, lo Bs = (0.50±0.05) GeV (for P° meson) have 
been used. It is clear that these hadronic parameters give the dominant theoretical uncertainties. 
They quantify the SU(3)-symmetry breaking effects in the form factors in pQCD a ppro ach. To 



make a comparison, we also collect the results using other approaches 0,1(3, 8, 39, 40]. Prom 
table llV] we can see that most of our results are consistent with others within theoretical errors. 
Likewise, B — > A form factors are defined by: 

2iA{q 2 ^ 



(A(P 2 ,e*)\q 1 ^ 5 b\B(P 1 )) 
(A(P 2 ,e*)\q^b\B(P 1 )) 

(A(P 2 ,e*)\qa^ 5 q v b\B(P 1 )) 
(A(P 2 ,e*)\qa^q,b\B(P 1 )) 



e^etP lp P 2a , 
m B - m A 

-2myVb(g 2 )^V - (m B - m A )^i(g 2 ) 
Q 



C-q 



+V 2 (q 2 ) 

m B - mA 
-2T x {q 2 )e^^P lp P 2(T , 



{Pi + P2) 



m 



A n 



iT 2 (q 2 ) [{m 2 B - m 2 A )e*» - (e* • ? )(Pi + P 2 )"] 



-iTz{ q 2 W ■ q) 



1 1 
m B — m A 



(Pi + P 2 



(56) 



with a factor — i different from B — > V and the factor m B + my (m B — my) is replaced by m B — m A 
(msA-rriA)- Similar with B — > V form factors, the relation 2rriAVo = (m B —rriA)Vi — {m B A-mA)V2 
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TABLE IV: B — ► V form factors at maximally recoil, i.e. q 2 = 0. The first error comes from decay constants 
of B mesons and shape parameters Lot] while the second one is from hard scale t and Aqcd- 



B 



B -> K* 



B 



B s -> if* 



LFQM[6] 



V 
A 
At 
A 2 



0.27 
0.28 
0.22 
0.20 



0.31 
0.31 
0.26 
0.24 



LCSR[8] 



V 
A 
At 
A 2 
T 2 



0.323 
0.303 
0.242 
0.221 
0.267 



0.411 
0.374 
0.292 
0.259 
0.333 



0.293 
0.281 
0.219 
0.198 
0.242 



0.311 
0.360 
0.233 
0.181 
0.260 



0.434 
0.474 
0.311 
0.234 
0.349 



LQCD[5] 



[39] 



V 
A 
Ax 
A 2 
Ti 



0.35 
0.30 
0.27 
0.26 



0.24 



SCET LCQM[40] 



V 



A Q 




A x 




A 2 




Tx=T 2 




n 




V 


0.21 


A, 


0.25 


Ax 


0.17 


A 2 


0.13 


Tx=T 2 


0.20 


T 3 


0.13 



0.298 
0.260 
0.227 
0.215 
0.260 
0.184 



0.339 
0.283 
0.248 
0.233 
0.290 
0.194 



0.275 
0.240 
0.209 
0.198 
0.239 
0.168 



0.323 
0.279 
0.228 
0.204 
0.271 
0.165 



0.329 
0.279 
0.232 
0.210 
0.276 
0.170 



This work 



-o 
f() 

-0 

fO 
-0 



05+0.00 
04-0.00 

05+0.00 
.05-0.01 

04+0.00 
.03-0.00 

03+0.00 
02-0.00 

04+0.00 
04-0.00 

03+0.00 
.02-0.00 



0.25^; 
0.301S 



0.19 
0.14 



+o. 



+0. 



0.23±g; 



0.13 



+o. 



05+0.00 
05-0.00 

06+0.00 
05-0.01 

04+0.00 
03-0.00 

03+0.00 

03- 0.00 

05+0.00 

04- 0.00 

03+0.00 
02-0.00 



0.20 
0.24 
0.15 
0.13 



+o 



+o 



0.1812 



0.12 



+o 



04+0. 
04-0 



05+0. 
04-0 



03+0. 
03-0 



03+0. 
02-0 



04+0. 
03-0 



03+0. 
02-0 



0.21 
0.26 
0.16 
0.12 
0.19 
0.11 



-0.04+0.00 
-0.03-0.01 



-0.05+0.00 
-0.04-0.01 



-0.03+0.00 
0.02-0.01 



-0.02+0.00 
-0.02-0.01 



-0.04+0.00 
-0.03-0.01 



-0.02+0.00 
0.02-0.01 



0.25^ 
0.30 
0.18 
0.13 



1.05+0.00 
1.04-0.01 



+ 



05+0.00 
.05-0.01 



+ 



03+0.00 
03-0.01 



02+0.00 
.02-0.01 



0.22±g 



04+0.00 
.04-0.01 



0.12 



+ 0. 



02+0.00 
02-0.01 



is obtained at q 2 = 0. In pQCD approach, B — > A form factors' formulas can be derived from the 
corresponding B — > V form factor formulas in eq. (|5QII55|) using the replacement in Eq. (139j) with 
the proper change of the momentum fraction. 

In the following, we will use a± to denote ai(1260), b\ to denote bi(1235). In Table W\ we give the 



17 



TABLE V: B — ► A form factors at maximally recoil, i.e. q 2 = 0. Results in the first line of each form 
factor are calculated using 9k = 45°, 9i Pl = 10° or 9s Pl — 38° , while the second line corresponds to the 
angle 9k = —45°, 9i Pl = 45° or 63 Pl = 50°. The errors are from: decay constants of meson and shape 
parameter w&; Gegenbauer moments in axial- vectors' LCD As. 



B -> ifi(1270) 



B -> hi(1170) 



S s -► JCi(1270) 



S -> /n(1380) 



/ii(1170) 



B s -> /ii(1380) 



.1 



-0.05 



+0.01+0.05 
-0.01-0.05 



0.33 



+0.07+0.05 
-0.06-0.05 



0.13] 
0.141° 



3.03+0.02 
3.02-0.02 

.03+0.02 
03-0.02 



-0.05 



+0.01+0.04 
-0.01-0.04 



0.07 



+0.01+0.01 
-0.01-0.01 



0.07 



+0.01+0.00 

-0.01-0.03 



n 1 -1 +0.02+0.04 
-U.11_q .02-0.04 



n no +0.00+0.00 
U.U/_q oi-O.OO 



-0.03 



+0.01+0.04 
-0.01-0.04 



-0.17 
-0.19 



03+0.02 
03-0.02 

03+0.02 
03-0.02 



Vq 



n -, 1 +0.02+1.10 
u - 11 -0.02-0.10 



0.60 



+0.12+0.10 
-0.11-0.10 



0.291°,; 
0.32±g; 



06+0.02 
05-0.02 



06+0.03 
06-0.03 



_n 1 O+0.02+0.09 
u - iz -0.02-0.09 



0.16 



+0.03+0.01 
-0.03-0.02 



0.11 



+0.03+0.09 
-0.02-0.09 



-0.05 



+0.01+0.00 
-0.01-0.00 



n -ic+o.03+0.01 

u - lu -0.03-0.01 

n ns+°- 01 + - 01 

u - uo -0.01-0.01 



-0.40±g 



-0.43 



+0 



06+0.03 
07-0.03 



07+0.03 
08-0.03 



-0.09 



+0.02+0.08 
-0.02-0.08 



n r 7 +0.12+0.08 
u - d ' -0.10-0.08 



0.22±g; 
0.24l° 



04+0.03 
04-0.03 



05+0.03 
04-0.03 



-0.09] 
-0.201° 



.01+0.07 
'.01-0.07 



04+0.07 
04-0.07 



0.12 



+0.03+0.02 
-0.02-0.02 



0.12 



+0.02+0.01 
-0.02-0.02 



-0.31 



+0 



n ru+°- 01 + - 01 

-u.ut_ oi_o .01 



u - ud -0.01-0.01 



-0.34j]° 



05+0.04 
06-0.04 



05+0.04 
06-0.04 



Vo 



n in +o.oi+o.o2 

-U.1U_ .02-0.02 



0.031!] 



0.11 



+0.02+0.02 
-0.02-0.02 



0.04 



+0. 



01+0.01 
01-0.01 



01+0.01 
01-0.01 



-0.01 
-0.17 



0.02 
0.02 



03+0.02 
04-0.02 



0.01 

-0.004 



ho. 00+0. 01 
-0.00-0.01 



+0.001+0.001 
-0.001-0.001 



n no+o.oo+o.oo 
u - uz -o.oo-o.oo 



-0.021°) 



-0.01 



+0.00+0.00 
-0.00-0.00 



-0.02 



+0 



00+0.01 
00-0.01 



00+0.01 
00-0.01 



Ti(T 2 ) 



n nfi+°- 01 + - 07 

u - uo -0.01-0.07 

n 4 r+0.10+0.07 
u -^°-0.08-0.07 



0.181°; 
0.201°,; 



04+0.02 
03-0.02 



04+0.02 
03-0.02 



-0.08l°) 



-0.12 



01+0.06 
01-0.06 

02+0.06 
03-0.06 



0.10 



+0.02+0.01 
-0.02-0.01 



n n/i+°- 06 + 04 

-U.W±_ Q oi-0 .00 



n in+o-02+o.oi 

u - lu -0.02-0.01 

_ n nA+0-00+0.01 
u - u ^-0.00-0.01 



-0.251° 
-0.27l° 



04+0.03 
05-0.03 



04+0.03 
05-0.03 



la 



n 1 1 +0.02+0.03 
-U.ll_ .03-0.03 



0.17 



+0.04+0.03 
-0.03-0.03 



0.05j 
0.06±g 



1.01+0.01 
1.01-0.01 



01+0.01 
01-0.01 



-0.02 
-0.19 



0.02 
0.02 



04+0.02 
05-0.02 



0.02 

-0.01 



-0.00+0.01 
-0.00-0.01 



+0.00+0.00 
-0.00-0.00 



n nq+0.01+0.01 
u - UJ -o.oi-o.oi 



-0.01 



+0.00+0.00 
-0.00-0.00 



-0.05^ 
-0.05 



'.01+0.02 
1.01-0.02 



+0 



01+0.02 
01-0.02 



B -> if 1 (1400) 



B -> /i(1285) 



B s -> i^i(1400) 



B-^ /1 (1420) 



B, 



A (1420) 



B, 



A (1285) 



A 



n n^+°- 01 + - 05 

u - ud -0.01-0.05 



0.34 



+0.07+0.05 
-0.06-0.05 



0.19; 

0.18 



04+0.02 
03-0.02 



04+0.02 
03-0.02 



n 1 O+0.03+0.04 
u - lz -0.02-0.04 

n n^+°- 01 + - 04 

U.UO_o.oi_o.04 



-0.011° 



-0.05 



-0 



00+0.00 
00-0.00 



01+0.01 
01-0.01 



-0.24;! 

-0.23 



i.04- 
1.04- 



0.03 
0.03 



-o.oil° 



+0 



04+0.03 
04-0.03 



-0.06 



+0 



00+0.00 
00-0.00 



01+0.00 
01-0.00 



Va 



_ n 19+0.02+0.11 

u - lz -0.03-0.11 

n R4+0.13+0.11 

u -°^-0.11-0.11 



0.261°, 



0.25 



05+0.06 
05-0.06 



05+0.06 
04-0.06 



n l7l +0.03+0.10 
u - Aft -0. 04-0. 10 

n 1 9+0.02+0.10 

u - iz -0.02-0.10 



-0.011° 

-0.071°, 



00+0.01 
00-0.01 



01+0.02 
01-0.02 



-0.341°,; 



-0.33 



+0 



05+0.08 
06-0.08 

05+0.08 
06-0.08 



-o.oil° 



-0.08 



00+0.00 
00-0.00 



01+0.01 
01-0.01 



Vs. 



0.09 
0.62 



02+0.09 
02-0.09 



13+0.09 
11-0.09 



0.32+° 



0.32 



+0. 



07+0.04 
06-0.04 



07+0.04 
06-0.04 



n 9/1+0.05+0.08 
u - zy: -0.04-0.08 

-0 nq+ ' 01 + - 08 

u - uy -0.01-0.08 



-0.021°, 



-0.09 



-0 



00+0.00 
00-0.01 



02+0.01 
02-0.01 



-0.44^ 
-0.43 



'.07+0.05 
1.08-0.05 



-0.021°, 



+0 



07+0.05 
08-0.05 



-0.10 



+0 



00+0.00 
00-0.00 



02+0.01 
02-0.01 



Vo 



0.101°; 
0.101°; 



02+0.02 
02-0.02 



02+0.02 
02-0.02 



O.IO; 
0.10; 



02+0.00 
02-0.00 



02+0.00 
02-0.00 



n 9n +o.o5+o.o2 

u - zu -0.04-0.02 

_n nn+ -o°+ 02 

u - uu -0.00-0.02 



-0.011° 
-0.021°, 



00+0.00 
00-0.00 



00+0.02 
01-0.02 



-0.11] 

-o.ir 



02+0.04 
02-0.04 



02+0.04 
02-0.04 



-0.011° 
-0.031°, 



00+0.00 
00-0.00 



00+0.00 
01-0.00 



Ti(T 2 ) 



-0. 



0.06 



0.481°, 



02+0.07 
01-0.07 



10+0.07 
09-0.07 



0.26 
0.25 



05+0.03 
05-0.03 



05+0.03 
05-0.03 



0.14 



-0.03+0.06 
-0.03-0.06 



n n o+o.oi+o.o6 

u - uo -0.01-0.06 



-0.01 
-0.07 



-0 



+0 



00+0.00 
00-0.00 



01+0.01 
02-0.01 



-0.33] 
-0.32 



'.05+0.04 
1.06-0.04 



-0.01 



+0 



+0 



05+0.04 
06-0.04 



-0.081°, 



00+0.00 
00-0.00 



01+0.01 
01-0.01 



la 



0.131°; 
0.151° 



03+0.03 
02-0.03 



03+0.03 
03-0.03 



0.14 
0.13 



03+0.04 
03-0.04 



03+0.04 
02-0.04 



n 9/1+0.06+0.02 
u - zq: -0.05-0.02 



-0.01 



+0.01+0.02 
-0.00-0.02 



-0.01 
-0.03 



+0 



-0 



00+0.00 
00-0.00 



01+0.00 
01-0.00 



-o.i6l°; 

-0.15 H 



03+0.00 
03-0.00 

'.02+0.00 
1.03-0.00 



-0.011° 
-0.041° 



00+0.00 
00-0.00 



01+0.00 
01-0.00 



This 



B —fbi 



B 



LFQM[6,_7] 



B^bx{K lB ) 



B -> clx{Kxa) 



A 

V 
Vx 
V 2 
Tx(T 2 ) 
T, 



0.191°; 

0.451° 

0.331° 
+0. 



0.03 
0.27 
0.06 



+0. 



+0. 



04+0.03 
03-0.02 

09+0.04 
.08-0.04 

07+0.04 
06-0.04 

01+0.01 
.00-0.01 

06+0.03 
05-0.03 

01+0.02 
01-0.02 



0.26 



+0. 



0.34l° 
0.43l° 
0.141° 
0.34l° 
0.191° 



06+0.03 

05- 0.03 

07+0.08 
.06-0.08 

09+0.05 
08-0.05 

03+0.00 
.03-0.00 

07+0.05 

06- 0.05 

04+0.01 
03-0.01 



A 

Vo 
Vx 
V 2 
Tx(T 2 ) 
T, 



0.10 (0.11) 
0.39(0.41) 
0.18(0.19) 
-0.03(-0.05) 
-(0.13) 
- (-0.07) 



0.25(0.26) 
0.13(0.14) 
0.37(0.39) 
0.18(0.17) 
-(0.11) 
-(0.19) 
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TABLE VI: Distinct contributions to form factor T\ from various distribution amplitudes. 





B^p 


B -> ax (1260) 


B -> 6x(1235) 




0.086 


0.142 


0.080 


4> a 


0.047 


0.086 


0.086 




0.063 


0.115 


0.100 


total 


0.196 


0.343 


0.266 



numerical results for B — > A form factors, in which we have used minus values for decays constants 
of 1 P\ mesons 3 . The errors are from: decay constants of Br s \ mesons and shape parameters cj&; 
Gegenbauer moments in axial- vectors' LCDAs. In the calculation, we use the mass of the two 
physical states i\i(1270), i^i(1400) as that of two spin states K\a, K\b for simplicity and similar 
for the branching ratios which are given in the following. That only involves a slight difference 
to the form factors. As the quark contents (to be precise the mixing angles) of the axial- vectors 
K\(f\, hi) have not been uniquely determined, we give two different kinds of results for form factors 
as in Ref. (37]]: the results in the first line are calculated using 9 k = 45° while the second line 
corresponds to the angle 9k = —45°. This is also done for the results involving the flavoi'-singlet 
and flavor-octet mesons: the results in the first line are calculated using 9ip 1 = 10°, #3 p l = 38°; 
while the second line corresponds to the angle 9ip 1 = 45°, 9zp 1 = 50°. 
A number of remarks on B — > A form factors are in order. 

1. Form factors are strongly dependent on mixing angles. Many of them even are different by 
order of magnitude, because the mixing angles describe directly the inner quark contents 
of the meson. The large difference of form factors surely will induce large differences in 
branching ratios, which we will see later. 

2. We give a comparison of the B — > p, B — > ai(1260) and B — > 6i (1235) form factors. Form 
factors Vo, V±, T\ for B — > A transition are larger than the corresponding B — > V ones. It 
seems that the form factor j[ B ^( a iM) [ s somewhat equal to or even smaller than V B ^ P . But 
actually that is artificial: as in Eq. (|5U|) , the pre- factor is + my while for B — > A form 
factor A, the factor becomes mp — rriA- We take T\ as an example to explain the reason 

3 Decay constants given in QCD sum rules [36l . [37| are both positive for two kinds of axial- vectors and we find that 
this will give negative values for B Pi form factors, like in [4lT |. For non-strange 1 P\ mesons, this minus sign 
will not give any differences as it can not be observed experimentally. But we should point out the minus sign will 
affect the mixing between K\a and K\b by changing the mixing angle 9 to —9. 
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for the large B — > A form factors. In table IVIl we give contributions from three kinds of 
LCDAs: 4> T , 4> v and (f> a . The contribution from (f> T is larger for B — > a\, than the other two 
transitions only because the axial-vector a\ decay constant is larger. Furthermore, larger 
axial vector meson mass implies larger contribution from twist-3 distribution amplitudes 4> v , 
(p a for both of T^ bl and T^ ai . 

3. In our calculation for form factors involving f\ mesons, we have used the mixing angle 
between the octet and singlet: 9 = 38° (50°) which is very close to the ideal mixing angle 
6 = 35.3°. That implies that the lighter meson /i(1285) is almost made up of uu +^ d while 
the heavier meson /i(1420) is dominated by the ss component. Thus B — > /i(1420) and 
B s — ► /i(1285) form factors are suppressed by the flavor structure and are numerically small. 
The form factors involving hi are similar if the mixing angle is taken as 45° . 

4. From the table [V] we can see that the form factor T^^ ai is almost equal to T± . In the 
flavor SU(3) symmetry limit, B — > K\a and B — > K\b form factors are also almost equal with 
each other. But the physical states ATi(1270) and A'i(1400) are mixtures of B — > K\a,\b- 
With the mixing angle Ok = ±45°, the B d)S — > Ai(1270)(Ai(1400)) form factors are either 
enhanced by a factor y/2 or highly suppressed. This feature will definitely affect branching 
ratios which will be discussed later. 

Up to now, there are not too many experimental constraints on B — > A form factors. However 



there are lots of studies using some non-perturbative met 



43, 



38 



45, 



rods: quark meson model [42|, ISGW 



4q ] and light-front quark model 



441 ], QCD sum rules and light-cone sum rules 
Results in LFQM are also given in table |V] to make a comparison. These two approaches are 
very different in the treatment of dynamics of transition form factors, but at first we will analyze 
differences caused by non-perturbative inputs. For B — > a\ and B — > K\a form factors, most of 
our results (except Vq and T\^) are slightly larger than or almost equal with these of evaluated in 
LFQM, as slightly larger decay constants for a\ and K\a are used: f ai = 203 MeV and fx 1A = 186 
MeV. The form factor Vq is calculated by the relation 2rriAVo = (tub — mA)V\ — {tub + rn,A)V2- 
Small differences in V\ and V2 have induced a large difference in Vo, which could be reduced in 
future studies using more precise hadronic inputs. We have found there are large differences in 
B Pi transition form factors. As the decay constant of b\ is zero in isospin limit, thus in LFQM 
the shape parameter u can not be directly determined and Cheng and Chua used the same value 



with that of a± 



6J. It is also similar for K\b'- they used the same shape parameter with that of K\a 



which predicts fx 1B = 11 MeV. Compared with the QCD sum rule results /k 1b = fx B x % given 
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in table HH and IIIIl we can see: although they are consistent within large theoretical errors, there are 
still large differences in the central value. Thus our predictions for B — ► l P\ form factors (central 
values) are larger than those in LFQM. We have to confess that the differences in decays constants 
are not responsible for all differences in form factors. That may arise from further differences in 
the dynamics. Compared with the recent light-cone sum rules results [461 ]: 



V^ Ol (0) = 0.3031 



0.035' 



v Bhl (o) 



-0.356 



+0.039 
-0.033' 



Vq Kia (0) 
V BKlB {0) 



316+ - 048 

u -° lu -0.042' 



-0.360 



+0.030 
-0.028' 



(57) 
(58) 



where uncertainties are from Borel window and input parameters, we can see that they are well con 



sistent with our ca 
of B° 



culations in pQCD approach. As mentioned in Ref. 



a t 7T [47J favors V ^ ai (0) ~ 0.30 and this is very close to our result: V 



411 ] , the Babar measurement 
B ^ ai (0) = 0.34. 

It is also noted that there are large uncertainties in our numerical results, especially due to large 
uncertainties of hadronic parameters, such as large Gegenbauer moments uncertainties shown in 
table EH 



IV. CALCULATION OF RADIATIVE DECAY B — > Vj 



A. The factorization formulae for decay amplitude 



For convenience, we define a common factor F which appears in many diagrams by: 

_ em 5 B C F 



7T 



(59) 



As we have mentioned in the above section, we have to use the amplitudes with distinct chiralities. 
The explicit factorization formulae for the left-handed and right-handed photon from operator O-j^ 
depicted in Fig. Q] are given by: 

rl poo f 

= 4r 6 F J dxidx 2 J b 1 db 1 b 2 db 2 (j)B(xi,bi)lc 7l (t a )E e (t a ) 



[1 + x 2 )(pv(x2) + rv(l - 2x 2 )(<Pv(x 2 ) + 4>v{xi)) h e (xix 2 ,x 2 ,b 1 ,b 



+ry 



(Pv(x 2 ) + (pv(x 2 ) C 7l (t' a )E' e (t' a )h e (xix 2 ,xi,b 2 ,bi 



n 1 



(60) 
(61) 



where the left-helicity amplitude is from the m& term in the effective Hamiltonian and the right- 
helicity amplitude is from the mo term which is obviously highly suppressed. This Of^ contribution 
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D b 



Q 1 



D 



7 7 

FIG. 3: Feynman diagrams of the chromo-magnetic penguin operator 0% g 



is the dominant one characterizing by the form factor Ti(T 2 ). The formulas in eq, (|6UI61|) are the 
same as that in eq. (|54p times the Wilson coefficient Cj. 

In pQCD approach, a hard gluon is required to kick the soft spectator in the B meson to 
turn into an energetic collinear anti-quark. This gluon could be generated from QCD interaction 
Hamiltonian or from the Og g operator. In Fig. O we give the four diagrams from Og g operator 
given in the effective Hamiltonian. The factorization formulae for the first two diagrams in Fig. [3] 
are 



M 



L(a) 



2r b F J dxidx 2 J bidb 1 b 2 db 2 (l)B{xi,b 1 )lQ D Csg{tb)E e (t b ) 

(2x 2 ~ xi)4>y(x 2 ) - 3x 2 r v (cj)y(x 2 ) + 4> v { x 2)) h e (x 1 x 2 ,x 2 - l,h,b 2 



+ 



xi^ v (x 2 ) + x 2 r v (^ v (x 2 ) + 0v(x 2 )) QbC 8g (t' b )E' e (t' b )h e (xiX2, 1 + x 1} b 2 , h) 



M 



R(a) 



(62) 



(63) 



If we consider the last two diagrams in Fig. O there will be more sources to generate the right- 
handed photon in addition to the mo term in the effective electro- weak Hamiltonian. The third 
diagram can give a small contribution which is from the higher twist component: 

Ms 9 W (Q 9 ) = 2r bQg F J o d Xl dx 2 J^ b 1 db 1 b 2 db2Mxuh)lc 8g (t c )E e (t c ) 

x (2 + x 2 - xi)4>y(x 2 ) + 3x 2 r v (<p v (x 2 ) + 4>v( x 2)) h e (x\ - x 2 , -x 2 , b±, b 2 ) 
+ x 2 r 2 ((j) v (x 2 ) + 4> v {x 2 )) - xi4>y(x 2 ) C 8g (t' c )E' e (t' c )h e (x 1 - x 2 ,x 1 ,b 2 ,b 1 ) 

rl f'OO 

-2r D Q q F I dx\dx 2 / bidbib 2 db 2 <j)B{xi,h)C$g(t c )E e (t c ) 
Jo Jo 

x 3x 2 r v ((j)y(x 2 ) - (j) v (x 2 )) h e (x\ - x 2 , -x 2 , &i, 62), (64) 
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D 



O, 



d b 



D 



FIG. 4: Feynman diagrams in which the operator O2 is inserted in the loop with a photon emitted from the 
external quark line 



M^ b \Q q ) 



rl poo r 

2r D Q q F J dxidx 2 J bidb 1 b 2 db 2 (l)B{xi,bi)lCsg{t c )E e (t c ) 
(2 + x 2 - xi)4>y{x 2 ) + 3x 2 r v ((j)y(x 2 ) + 4>v{x2)) h e (xi - x 2 , -x 2 , h,b 2 ) 



+ 



x 2 r v ((p v (x 2 ) +(f> v (x 2 )) - xi(f> v (x 2 ) C 8g (t c )E e (t c )h e (xi - x 2 ,xi,b 2 ,b 



+2r b Q q F / dx\dx 2 I bidb]b 2 db 2 (t)B{xi,bi)C%g{t c )E e (t c ) 
Jo Jo 

x 3x 2 r v ((fiy(x 2 ) - (j) v (x 2 )) h e (xi - x 2 , -x 2 ,bi,b 2 ). 



(65) 



Next we want to mention some higher order corrections as usual which may give important 
contributions: charm and up quark loop (0(a s )) contributions in Fig. 01 Fig. [5] and Fig. [6l It 
should be pointed out that these contributions are not related to next-to-leading order corrections 
in pQCD approach, while next-to-leading order corrections to the exclusive processes 7r7* — > 7 in 
pQCD approach have been investigated in Ref. [l^ . 

We use the subtitle "quark line photon emission" to denote that a photon is emitted through 
the external quark lines as in FigJU We define the c and u loop function in order that the b — > Dg 
vertex can be expressed as 1)7^(1 — 7 5 )/^ A ai/ b. It has the gauge invariant form [48j as follows: 



ja 



T^-ji^g^u - fe/iM / dxx(l - x) 



8tt 2 



G{mlk\t) 



1 + log 
2 



m? — x(l — x) 



(66) 



where k is the gluon momentum and rrii is the loop internal quark mass. G is the function from 
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the loop integration: 



G(mf,k 2 ,t) 



5 4mf , to~ 

3 + ^- ln 7^ 



2mi 



+e{k 2 )e{Am 2 l -k 2 )- 



k 2 

5 Am 2 , to, 2 
3 + ^- ln # 



4m 2 j i/l - 4m 2 /k 2 - 1 



/v 2 



^1 - 4TO 2 /fc 2 + 1 



-2 1 



2m\ 
~~k? 



+9(k 2 -4m 2 )- 



2m' 



Am 2 

„ — 1 arctan 



y/4m?/k 2 - 1 



5 4771 2 , TO 2 



A' 2 



4m| 
fc 2 



In 



1 - v /l-4TO 2 /fc 2 



(67) 



1 + y/l - 4TO 2 /fc 2 

The loop function G has the dependence of gluon momentum square of k 2 . But there is no 
singularity when we take the limit of k — ► 0, so we can neglect kx components of k 2 in the loop 
function G. Using this effective vertex, the factorization formulae for the first two diagrams of 
Fig|4]is calculated as: 



M 



L(a) 
li 







Q D F / dx 1 dx 2 / b 1 db 1 b 2 db 2 <j)B(xi,bi)C2(tb)E e (t b ) G(m 2 , -xix 2 m 2 B ,t b ) 



x 2 r v ((t> v {x 2 ) +<j)y{x 2 )) + 3x 1 x 2 cf> v (x 2 ) h e (xix 2 ,x 2 - l,b 1 ,b 2 ), 
M R{a) = _Q h p j dxidx 2 j bidbib2db 2 (f>B(xi,bi)C2(t' b )E' e {t! b ) G(m 2 , -xiX 2 m B , t' b ) 
xx 1 x 2 r v ((j)y(x 2 ) - <p^(x2))h e (xix 2 , 1 + xi,b 2 , &i), 
where Qd = — , . For the other two diagrams of FigJH the factorization formulas are 



(68) 



(69) 



M L J b) (Q q ) 



rl roo ( 

Q q F j dx\dx 2 j bidbib 2 db 2 4> B {xi,bi)lC 2 (t c )E e (t c ) 
x 2 r v {\ + 2x 2 )[^(x 2 ) + 4fy{x 2 )\ + 3(x 2 - xi)^(x 2 ) 



G(m 2 ,-(xx -x 2 )m 2 B ,t c ) - - h e {x\ - x 2 , -x 2 , b\, b 2 ) 



+ 



Mf} b) (Q q ) 



x 2 r v {4>v{x 2 ) + <Pv(x 2 )) - xi(py(x 2 ) G(m 2 , -{x\ - x 2 )m 2 B ,t' c ) 
x C 2 (t' c )E' e (lf c )h e (xt - x 2 , xx , 6 2 , h) j , 

Q q F dxidx 2 bidbib 2 db 2 4>B{xiM)C 2 {t c )E e {t c )h e (x 1 - x 2 , -x 2 ,bi,b 2 
Jo Jo 

2i r 
G(m 2 , -(xi - x 2 )m 2 B ,t c ) - - x (2 + x 2 )x 2 r v <ffy (x 2 ) - (j)y (x 2 ) 



(70) 



(71) 



In Fig. [5l we give the diagrams in which a photon emitted from the internal loop quark line. 



The sum of the effective vertex b — > D^g* in FigJS] has been derived by [4 

/ = DY(l - j 5 )T a bI^ p A^A a \ 



(72) 
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FIG. 5: Fcynman diagrams in which the operator 2 is inserted in the loop with both of the photon and 
the virtual gluon emitted from the internal quark line. 

with the tensor structure given by 



= At [(q ■ k)e^p a (q - k) a + e upaT q a k T kp - e fipar q a k T q u ] 
~l~-<4-5 l^nparQ k k v — k € pV p a q ] , 



(73) 



and 



A 5 



Meg f 1 ^ I 1 x rf xy 

3vr 2 Jq Jq x(1 — x)k 2 + 2xyq ■ k — mf + is ' 

Iq Jq x(1 — x)k z + 2xyq ■ k — mf + ie 



3tt 2 



(74) 
(75) 



where q is the momentum of the photon q = P B — Pv, and k is the momentum of the gluon. Then 
the amplitudes can be expressed as follows: 



Mi 



> rl rl—X rl 

■F I dx I dy dxidx 2 \ b\db\$ B (xxM)C 2 {t d )a s (t d ) exp[-S B (t d )] 
' Jo Jo Jo 



K 

2 2 x ^ x y x z 

xyx 2 m B — mi 



(1 - x 2 )r v (<t> v v (x 2 ) + <^(x 2 )) - (1 + 2x 1 )cj ) 1 v (x 2 ) 



+x(l - x) xlrv((pv(x 2 ) + 4>v(xi)) + ^XlX 2 (j)y{x 2 ) 



-F J dx J dy J dx±dx 2 \ bidbi^ B (xi,bi)C 2 (t d )a s (t d ) exp[-S B (t d )] 



K 2 

: 2 2 X X V X 2 r V 

xyx 2 m B — m\ 



where the function h' e is defined by: 



h' e = K Q y^xim B bx)- [^(5 2 )iiCo(6 1 v1^|)+^(-B 2 )^ J ffo(&iv^|) 



with 



(76) 



(77) 



(78) 



B 2 = x\x 2 m\ 



y 



-x 2 m B + 



(79) 



1 — x D x(l — x) 

Diagrams in Fig. [6] in which the photon is emitted from the external loop and the gluon is 
emitted from QCD interaction Hamiltonian do not give any contribution to B — > V7. The 
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FIG. 6: Feynman diagrams in which the operator 2 is inserted in the loop with a photon emitted from the 
internal quark line 




7 K 

FIG. 7: Feynman diagrams for annihilation topologies 



reason is as follows. Similar with b — > Dg, the vertex for b — ► D7 can also be expressed as 
5/^ — kpki^b only with a different coefficient. For an on-shell photon, the fol- 
lowing conditions are required: k 2 = and e - k = 0, thus the contribution from diagrams in Fig. [6] 
vanishes in b — > s(d) r y decays. But it should be noted that these diagrams can give a non-zero 
contribution to b — > s(d)j* — > s{d)l + l~ . 

In annihilation diagrams, there are three different kinds of operators in the (8) depicted in Fig. [71 
In the following, we use LL to denote the left-handed current between b and q quark and the left- 
handed current between the final state two quarks; LR denotes the left-handed current between b 
and q quark and the right-handed current between the final state two quarks; we use SP to denote 
the (S — P)(S + P) current which is from the Fierz transformation of (V — A){V + A) operators. 
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The factorization formulae for these diagrams are given by: 



•^■ann (fliiQqi) ~ -^ann \ a iiQ 



<n . 



Jvl ann \ a t^ 



31 - 



"B JO 
^■ann ( a itQqi) 

F 3V6Q b r v f V 7T f 1 ^ 



(81) 



2m| 



bidb 1 a i (t e )E a (t e )4>B(xi,b 1 )K (y/l + xim B bi 



■^■ann ' '( a «i Q 92' Q 93) ~~ -^oin ^ { a i 1 Q q2 i Q qi, 



2m| jo 



^(x 2 ) + ^(x 2 ) 



x <|Q g2 a J (t/X(i/)i|^ 1) (v / r 3 ^m i? 6 2 ) 



-x 2 Q qi ai(t' f )E' a (t' f )i^H^\yfx^m B b 2 ) 



■M a nn \ a i 1 Q '132 ' Q 93) ~~ -^onn ( a i ' Q<22 > Q<23 ) 



2m| jo 



x (I-12) 



(82) 



0v(^ 2 ) - ^(x 2 )]Q (?2 a i (t / )K(i/)^^ 1) (v / r^"ii?fe2) 



<^(x 2 ) - <fr(x 2 ) Q q Mt'f)K{t'f)i\H£X^x- 2 m B b 2 ) 



(83) 



m 



B JO 



<£c 2 / b 2 db 2 cpy(x 2 ) 



(84) 



Q ?2 di(i/)K(*/)^- ff o 1) (Vl - ximBbi) + Q g3 a i (t / / )-E^(t / / )i|^ 1) ( % /i^m B 6 2 )|. 

Finally, there is another kind of contribution from : the neutral vector meson qq is generated 
by a photon as depicted in Fig. [8l Although these diagrams are suppressed by the electromagnetic 
coupling constant, the enhancement factor m#/ 'A.qcd can make it important in some cases 510 . 
We include these diagrams in our calculation. The first two diagrams of Fig. [8] are equal to each 
other, so are the last two diagrams. The factorization formulae are given by: 

3y/6a em Q q Q b r b f v f l 



M L en {Q q ) 



niBmv 



dxi I bidb\4>B{x\M) 
o Jo 



M*(Q q ) = ~M^{Q q ). 



x i C 7l {t e )E a (t e )K {^/TTx^m B bi) + C 7l (t' e )E a (t' e )Ko{^m B bi) \, (85) 



(86) 



r D \ aL 



27 




FIG. 8: Fcynman diagrams with double-photon contributions 
B. Numerical results of Branching ratios 



With those decay amplitude formulas for different Feynman diagrams in the last subsection, it 
is easy to get the total decay amplitude for each channel of B — > V7: B~ — > p~7, -B — ► p°7, 
5° cj 7 , B~ -> if*-7, £° ^7, B° -> A*° 7 , B° -> A'*° 7 , B° p° 7 , B° ury, and £° <fyy. 
The explicit expressions are shown in Appendix [Bj The CP-averaged decay width is then 

where the summation on polarizations is implemented. 

For CKM matrix elements, we use the same values as in Ref. \)\ '- 

\V ud \ = 0.974, \V US \ = 0.226, \V uh \ = (3.68±g;J|) x KT 3 , 

\Vtd\ = (8.20lg;p) x KT 3 , = 40.96 x 10~ 3 , \V tb \ = 1.0, (88) 
« = (99l£ 4 )°, 7 = (59.0l|:?)°, arg[-V^] = 1.0°, 

□ 

where we have adopted the updated results from [52| and drop the (small) errors on V u d, V us , Vt s 
and Vtb- The CKM factors mostly give an overall factor to branching ratios. However, the CKM 
angles do give large uncertainties to branching ratios of some decay modes and to all the non-zero 
CP asymmetries which will be discussed in the following subsection. 

CP-averaged branching ratios of B — > Vj decays are listed in table IVlTl The first error in these 
entries arises from the input hadronic parameters, which is dominated by B(B s )-meson decay 
constants (taken as f B = (0.19 ± 0.02) GeV and f Bs = (0.23 ± 0.02) GeV) and B (B s ) meson 
wave function shape parameters (taken as ujb = (0.40 ± 0.05) GeV and ub s = (0.50 ± 0.05) GeV). 
The second error is from the hard scale t, defined in Eqs. (Al) - (A8) in Appendix A, which 
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TABLE VII: CP-averaged branching ratios (xlO -6 ) of B — * V7 decays obtained in pQCD approach (This 
work); the errors for these entries correspond to uncertainties in input hadronic quantities, from the scale- 
dependence, and CKM matrix elements, respectively. For comparison, we also listed the current experimental 



measurements 
and in Ref. 6 Z 



I l| and theoretical estimates of branching ratios recently given in Ref. [3l| (QCDF) 
(SCET). 



iviocics 


or" fit? 




This work 


Exp. 


D — > l\ 7 


OO. O ± lO.O ± O.O 


/l£i_l_19_l_/l_l_9_l_1 


or 0+17.6+5.4+1.1 

JJ -°-12. 8-4. 0-1.1 


40.3 ±2.6 [HFAG] 




54.2 ± 13.2 ± 6.7 


43±11±4±2±1 


00 1 +17.3+5.5+1.1 
°°- 1 -12. 7-3. 8-1.1 


40.1 ± 2.0 [HFAG] 


B°s - 07 


39.4 ± 10.7 ±5.3 


43 ± 11±3±3± 1 


or n + 13.7+4.9 + 1.1 
od -°-10. 3-3. 5-1.1 


STtStg [Belle] 


Modes 


QCDF 


This work 


Exp. 


B~ -> p-j 


1.16 ±0.22 ±0.13 


1 1 r+0. 57+0. 18+0. 17 
x - ±d -0.39-0.11-0.09 


1.10±g;|| ± 0.09 [BaBar] 


0.55t°;™ [Belle] 


B° -> p° 7 


0.55 ±0.11 ±0.07 


n c.7+0. 26+0. 09+0. 08 
u -° ' -0.19-0.06-0.04 


0.791:^20 ± °- 06 [BaBar] 


1-25±8S1S5S [Bdlc] 


B° -> Lj 7 


0.44 ±0.09 ±0.05 


n ci +0.23+0.08+0.08 
u - dl -0. 17-0. 05-0. 03 


- 40 -o:lo ± °- 05 [BaBar] 


0-56t°; 2 ™o' [Belle] 


5° -> K*°-y 


1.26 ±0.25 ±0.18 


1 11+0.42+0.15+0.16 
iiil -0. 32-0. 12-0. 07 




B° -> (fry 




In r+2.8+2.1+l.l\ iq-6 
{' -°-2.1-0.9-0.5^ x iu 


<0.85 [HFAG] 


B° s -> p°7 




( I 7+0. 4+0.1+0. On 1n -3 
I 1 - '-0.4-0. 1-0. lJ x iu 




B° -> w 7 




fi o+0.4+0.1+0.1n i n -4 
U-^-O^-O^-O.lJ x iu 





we vary from 0.75t to 1.25t (not changing 1/bi), and from A.q^, d = 0.25 ± 0.05 GeV. This scale- 
dependence characterize the size of next-to-leading order contributions in pQCD approach. A part 
of this perturbative improvement coming from next-to-leading order Wilson coefficients is already 
available [22|. However, the complete next-to-leading order corrections to hard spectator kernels 
are still missing. The third error is the combined uncertainties in CKM matrix elements and angles 
of the unitarity triangle. It is clear that the largest uncertainty here is the first one from the input 
hadronic parameters. 

These ten B — > Vj decay channels can be divided into three different types: b — > s transitions, 
b — ► d transitions and purely annihilation decays. The first type contains B° — ► K*°j, B~ — > K*~j 
and B s — > (fry. Among these decays, the dominant contribution from Oj^, is proportional to 
VtbV^ ~ A 2 . This contribution can be related to the form factor T^^ v . In the flavor SU(3) 
symmetry limit, form factors for these three channels should be equal which could also relates the 



three decays. We do obtain similar branching ratios for this kind decays only with small deviations. 



The penguin contribution in b — > d'y processes is proportional to VtbV.*^ ~ A 3 , which is expected 
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to be suppressed by one order magnitude relative to the b — > s transitions. In our calculation, 
branching ratio results for pure annihilation processes are mainly from two-photon diagrams in 
Fig. Thus B1Z(B — ► cfyy) is surely consistent with Ref. [5l|. This feature can also certainly 
interpret the large differences between B s — ► p°j and B s — > W7. This contribution is proportional 
to the charge of the constitute quark. This factor is | — ^ for p while § + -j" for w. Thus the 
branching ratio of i? s — ► p°j is one order in magnitude larger than that of (B.^ -7 toy)- 
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nt 



59|, 



re literature, there are many studies concentrating onB^ V7 [la, [l7| 



18 



53 



54 



55 



56 



60, 



61 



62 



63 



57|, 



64 ] . Recently, a comprehensive study [31]] using QCDF method and QCD sum 



rules appears. In that paper, the authors used QCDF approach to calculate all B — ► Vj approach 
and included some power corrections: weak annihilation contributions, the soft-gluon emission from 
quark (charm and light-quark) loops and long distance photon emission from the soft quarks. In 
Table IVUl we quote them to make a comparison. Their uncertainties come from form factors, the 
renormalization scale, the soft-gluon terms, the CKM parameters, decay constants, Gegenbauer 
moments, the first inverse moments of B(B S ) mesons and quark masses, etc. Their results for 
branching ratios of B~~ — > K*~^ and B° — ► K*°j are about (20—50)% larger than our predictions. 
The main reason is differences in form factors T\\ they used Tf^ K = 0.31 ± 0.04 while our 
calculation gives a smaller T x = 0.23±g;g|. The smaller form factor is also preferred by recent 
Lattice QCD result: Tf~* K * = 0.24 ± 0.03 j^. Although the difference is not too large, it can 
already induce a sizable difference to branching ratios, b — > d transitions and B s — ► (jry are well 



311 ] are smaller and almost 



consistent with each other, as the effective form factors used in Ref. 
equal to our results. Very recently, the authors in Ref. [64| used the more theoretical approach 
SCET to investigate the three b — > s decays channels: B~~ — ► K*~j, B° — > K*°j and B s — ► (jyy. 
After integrating out the hard scale which results in SCET/, contributions to B — > V7 decay 
amplitudes can be divided into two different groups: contributions from operators J a and Jb- 
The .B-type operator is power suppressed in SCET/ but it can give leading power contributions 
when matching onto SCET// as A-type operator receives power suppressions. When performing 
the matching from SCET/ to SCET//, we have to be cautious about the ^4-type operator as this 
term suffers from the end-point singularities. Thus one has to leave it as a non-perturbative free 
parameter determined from experiments or some non-perturbative QCD approaches, but recent 
studies using zero-bin subtractions show that this term can also factorized in rapidity space [65J] . In 
Ref. 64J], the authors calculated the two- loop corrections (a^) to short-distance coefficients of the 
A-type operator (called vertex functions) determined from QCD to SCET/ matching and utilized 
the physical B — > V form factor T\ to extract the soft form factor in SCET. While for the B- 
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type operator which can been factorized into convolutions of LCDAs and hard kernels, both of 
the Wilson coefficient and jet function have been calculated up to one loop order. Short-distance 
Wilson coefficients for the -B-type operator are extracted at the scale m& and have been evolved 
down the intermediate scale \i ~ 1.5 GeV using renormalization group equations. With these 
results, the authors find: compared with the leading order contribution, the order a s corrections 
from vertex corrections and hard spectator scattering can be as large as 10% and both of them 
also provide imaginary amplitudes about 5% in magnitude; the order a 2 corrections are not too 
large. We quote their final results for branching ratios in table IVIII and they are also consistent 
with ours. 

The three experimental collaborations, BaBar 



66j, Belle 



67fl and CLEO 



, have reported their 



measurements on B1Z{B — > K*j). Since all of these results are well consistent with each other, 
we quote the averaging results from Heavy Flavor Averaging Group (HFAG) |J in table rvTTl We 



691 ]. All of them agree with our calculations. 



also include the very recent result on BTZ(B® — ► <jyy) 
On the b — ► d transition, branching ratios of some channels have 

and Belle 7l( collaboration. We find that our results agree well with BaBar's results but not 
with Belle's central value results. In flavor SU(3) symmetry limit, the relation B1Z(B~ — > p~j) = 
2BK(B° -> p°7) = BK(B° -> wy) should be held. The small deviation in our calculation is caused 
by the SU(3) symmetry breaking effect. Since the electro-magnetic penguin operator Oj^ gives 
the dominant contribution, it is difficult to understand the results from Belle collaboration: why 
is the branching ratio of B° — > p°7 larger than B~ — > p~ 7. But before we conclude it is the signal 
for non-standard model scenarios, it is necessary to re-examine this channel on the experimental 
side. All other decay modes, including B s decays and annihilation type decays, have not been 
experimentally measured. 



C. CP asymmetry studies 

The direct CP-asymmetry in B — ► V7 is defined by: 

dir = BK(B -> Vj) - BU(B Vj) = lAjB^V^-lAiB^V^l 2 

CP BK(B ^Vj)+BK(B ^V~f) \A{B -> y 7 )| 2 + \A(B -> y 7 )| 2 ' 1 ; 

In order to give a non-zero direct CP asymmetry, we need two kinds of contributions with different 
strong phases and different weak phases. The magnitude of the CP asymmetry also depends on 
relative sizes of the two different amplitudes: if one amplitude is much larger than the other one, 
we can only get a small CP asymmetry. Since there is only penguin contribution in B° — > (fry 
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A d S>{B~ 


P 7) 


1 9 0+0.8+2.9+0.8 

— lz -°_o.3-l. 8-0.8' 


A%£(B° ■ 


-A) 


_ 19 ,,+0.2+1.8+0.5 
— -L^-4_o.4-2.4-0.9' 


A d &(B- - 




= -0.4 ±0.0 ±0.1 ±0.0 


A$£(B° -> 




n q+0.0+0. 0+0.0 
— u -°-0.0-0. 0-0.0' 


A^(B° 




19 1+0.0+1.8+0.5 
— 1Z.1_ 2 _ 2 4_ g, 


A&(B° a ■ 


-A) 


n 1+0.0+0.3+0.0 
— u - i -0.0-0. 1-0.0' 


A&(B° a - 


ir° 7 ) 


1 9 7+0.1+1.6+0.5 
— -0.5-2.3-0.9' 




-> UJj) 


n q+0.0+0. 9+0.0 
— u -°-0.0-0. 5-0.0' 


^(4° 


07) 


n q+0.0+0.0+0.0 

— u '°-0.0-0.0-0.0' 



process, the direct CP asymmetry is zero. In other B — > decays, there are contributions from 
the penguin operator and two kinds of tree operators (proportional to V c bV* dcs and KibV^ us )- 
Taking these amplitudes into account, we obtain the numerical results for direct CP asymmetries 
(in %) in other B — ► V7 decays as: 

(90) 
(91) 
(92) 
(93) 
(94) 
(95) 
(96) 
(97) 
(98) 

where the three kinds of errors are given as that of the branching ratios case. It is easy to see that 

theoretical uncertainties here are much smaller than branching ratios in table IVI!} especially the 

first one from hadronic input parameters, since they are mostly canceled in eq. (|89p . In the three 

b — > s channels B~ — ► K*~-y, B° — > K* ^/ and B s — > (^7, CKM matrix element for the magnetic 

penguin operator is VtbV* s ~ A 2 , while the tree operator is either proportional to V c b V* s ~ A 2 

° r KifeKTs ~ -^ 4 - The CKM matrix element in the first tree operator is almost parallel to the 

penguin operator. This kind of contribution has a same weak phase with the penguin contribution. 

The second tree operator is small in magnitude. Thus we expect small CP asymmetries in these 

channels. Experimentally, both BaBar and Belle collaboration give their combined measurements 

□ 

of the two B — > if* 7 channels [721 ]: 

f -1.5 ± 4.4 ± 1.2, [Belle], 
A$£(B^K*j) = { 1 J (99) 

[ -1.3 ± 3.6 ± 1.0, [BaBar]. 

Although the central value of CP asymmetry is larger than the one in our calculation in eq. (|92|93p , 
it is still consistent with zero. 

In annihilation- type decays B s — > p°(u)^, the tree amplitude is also suppressed by CKM matrix 
elements, thus the CP asymmetry is small too. In b — > d transitions B~ — > p*~7, B° — > p°(u>)j 
and B s — > K* "/, the CKM matrix element for the magnetic penguin operator is V^V^ ~ A 3 , while 
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the tree operator is either proportional to V c bV* d ~ A 3 or V u bV* d ~ A 3 . Then tree contribution is 
not suppressed and can be comparative with the penguin contribution. Thus we expect relatively 
large CP asymmetries in these four processes, which are also shown in the above. 

Restricting the final vector state V to have definite CP-parity, the time-dependent decay width 
for the B° — > f decay is: 

r(S°(f)->/) = e - r *r(5^/)[cosh(^)+iI / sinh(^) 

cos(Amt) - S f sin(Ami)] , (100) 

where Am = mu — tul > 0, T is the average decay width, and Ar = Th — is the difference 
of decay widths for the heavier and lighter B° mass eigenstates. The time dependent decay width 
T(B°(t) — > /) is obtained from the above expression by flipping the signs of the cos(Ami) and 
sin(Ami) terms. In the B d system, Ar is small and can be neglected. In the B s system, we expect 
a much larger decay width difference (AT/T)b = —0.127 ± 0.024 [3] within the standard model, 
while experimentally (AT/T)b 3 = — 0.33^q'5i [3], so that both Sf and Hf, can be extracted from 



the time dependent decays of B s mesons. The definition of the various quantities in the above 
equation follows: 



21m (1(A* l Al + A* r A R ) 



S f^ = ^WWW' (101) 



2 Re [UA* l Al+A* r Ar) 

Wi> - + < 102 » 

where A and A denote the amplitudes for the B and B meson decays, q/p is given in terms of the 
Bq-B® mixing matrix M12, 



M* 2 



M 



12 



(103) 



with 

<P d = - W g[{V t * d V tb f] = -2/3 , <t> s = - & Tg[{V t * s V tb ) 2 \ = 2e. (104) 

where the convention arg[V c fc] = arg[V cs ] = is adopted. 

In b — > Dj(D = d, s) processes, the dominant contribution to decay amplitudes comes from 
the chiral-odd dipole operator O7. As only left-handed quarks participate in the weak interaction, 
an effective operator of this type necessitates, a helicity flip on one of the external quark lines, 
which results in a factor mj (and a left-handed photon) in — > -Dl7l and a factor win (and a 
right-handed photon) in bi — > Dr^r- Hence, the emission of right-handed photon is suppressed 
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TABLE VIII: Mixing-induced CP-asymmctry parameters (in percentage) of B — > Vj decays obtained in 
pQCD approach. The errors are the same with table IVlIl The Hf parameter in B^ decays could hardly be 
measured as the decay width difference is small. 



Modes 


Sf 




B° -> K*°-y 


A n+0. 1+0.3+0.1 
^•"-0.0-0.2-0.1 


a i +0.0+0.2+0.1 
^-0. 1-0. 3-0.1 


B° a -> 07 


n Q+0.0+0.0+0.0 
"• z -0. 0-0. 0-0.0 


r r+0. 1+0.4+0.0 
°-°-Q. 2-0. 4-0.0 


B° -> p° 7 


n o+0. 1+0. 0+0.0 
"■°-0. 1-0.2-0.1 


n a+0, 3+0.5+0.1 
"■^-0.1-0.3-0.1 


B° -> w 7 


n .,+0.0+0. 2+0.0 
"■^-0.1-0.2-0.0 


n r +0. 1+0.4+0.1 
u - d -0. 0-0. 3-0.0 


B° s -> K* ^ 


n 7 +0.0+0.4+0.1 
"■ ' -0.1-0.2-0.1 


n O+0. 1+0.3+0.1 
u -°-0. 2-0. 3-0.0 



by a factor mo/mft. In the b — > D7 process, the emitted photon is predominantly left-handed, 
and right-handed in b decays. This leads to very small predictions of Sf and Hf. The mixing- 
induced CP asymmetry variables are calculated and summarized in table IVIIIi B° — ► K* ^ has 
been treated as an effective flavor eigenstate. Apparently, the numerical results agree with our 
expectations. On the experimental side, the mixing-induced CP asymmetries have been measured 
by Belle and BaBar as follows Q,H, 2S]: 



-0.79j#S ± 0.10, [Belle], 
S f (B^K*j^Ks^j) = \ (105) 

-0.08 ±0.31 ±0.05, [BaBar], 
S f (B° -> p°j) = -0.83 ±0.65 ±0.18. (106) 

They are consistent with zero since there are large uncertainties in these results. The theoretical 
results agree with the experimental data taking the experimental uncertainty into account. But as 
this parameter could be a good probe to detect the non-standard scenarios, more studies, including 
both of the precise experimental studies and the theoretical studies, are strongly deserving. 

D. Isospin asymmetry and U-spin asymmetry 

Apart from branching ratios and CP asymmetry variables, we will also consider some ratios 
of branching fractions defined below. In the evaluations for branching fractions, there are many 
uncertainties, especially from hadronic input parameters, which can blur our predictions, but we can 
improve the accuracy of our predictions by using ratios of branching fractions. Many uncertainties, 
such as those from decay constants, will cancel in these parameters. The most important ratios 
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are the parameters characterizing isospin asymmetries which are defined by: 

where the partial decay rates are CP-averaged. 

In the flavor SU(3) symmetry limit and if we neglect diagrams which are proportional to the 
quark charge, all of these three parameters should be equal to 0. The u meson decay constant 
is smaller than that for p° meson, the uu component contributes with a different sign and the 
electomagnetic diagrams with two photons are different in the charge factor. These differences 
make A(p, uS) deviate from (smaller than 0). The origins for deviations for Aj(p) and Aj(K*) 
from are similar: the spectator quarks are different and annihilation diagrams are also different. 
Taking on all those power suppressed contributions, our predictions are 

A( nuJ \ - _ n ii+o.oi+o.oi+o.oo (ua) 

A \P,U) — U.l±_ QQ_ 00 _Q QQ, ^IIUJ 

Mp) = o.oelHIlHn^ (in) 
MK*) = o.otf™£™£™. (112) 

As we expected, theoretical uncertainties due to the hadronic parameters are indeed smaller due to 
cancelations. Using experimental results listed in table IVTl] we can calculate the isospin asymmetry 
parameters from experiments 

' (». !') •!;;:!,; [BaBar], 
A(p,u) = { (113) 

-°- 55 -a27 [ Bel H 

-0.54+°-^ [BaBar], 
Mp) = { -°- 67 (114) 

3.89±H° 4 [Belle], 
Aj(K*) = 0.03 ±0.04, (115) 

where we have assumed all the uncertainties are not correlated by adding them quadratically. 
From the experimental results, except Aj(K*), we find that there are large differences between the 
results from the two collaborations. Our results are consistent with them, since the error bars in 
the experiments are too large. We wish a more precise measurement on these parameters in the 
future. 
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Apart from the isospin symmetry, U-spin symmetry is another kind symmetry which is well 
held in strong interactions. U-spin can connect two different kinds of weak decays 77, Q, Q] : 
b — ► s (AS = 1) and b — > d by exchange of d <-> s. The decay amplitudes of 6 — ► s process can be 
expressed as: 

A(B - /) = + WsA, (116) 

^(5 _> /) = KfeKTsX + K*04c, (117) 

while the decay amplitudes of UB — > £7/ are 

^(c/s _> [//) = y^Kdf/X + K* b u cd lU c , (H8) 

X(17S -> C//) = V ub V: d UA u + V cb V; d UA c . (119) 

Using the relation A u = UA U and A c = UA C in U-spin symmetry limit and the CKM unitarity 
relation 

lm(V: b V us V cd V c * s ) = -Im(V: b V ud V cd V c * d ), (120) 

we obtain 

\A(B - /)| 2 - - /)| 2 = -|^(17B - Uf)\ 2 + |^(Z7B - C//)| 2 . (121) 

This equation relates the differences of partial decay widths. The following radiative B decays 
can be related to each other by this symmetry: B~ — ► p~7 and B~ — > K*~^; B° — ► K*°~f and 
Bg — ► K*°j. As an example, we define the following parameter to test U-spin symmetry breaking: 

A = Acp(B~ - JT- 7 ) - Acf(S- - p" 7 ) x z^Sk^I) ■ (122) 
In our pQCD calculation, we find 

A = (-8.4t°;^±°;3) x 1(T 3 . (123) 

This result is close to 0. The U spin symmetry seems quite good here. But the most important 
reason is the small CP asymmetry in B~ — > K*~j and the small ratio ' absolute 

value for the b — > s channel's CP asymmetry is small. This direct CP asymmetry may be dramati- 
cally enhanced by some new physics with a different weak phase. The precise measurement of CP 
asymmetries can at least give a constraint on the non-standard model scenario parameters. 
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V. CALCULATION OF B — ► A(l 3 Pi) 7 AND ^(1^1)7 DECAYS 

The factorization formulae for B — ► j4 7 is more complicated than — > V'j because of the 
mixing between different mesons: K\a and K\b] fi and f$; h\ and /ig. The real physical states 
-fTi(1270) and i^i(1400) are mixtures of the K\a and K\b states with the mixing angle 9k- 

|#i(1270)) = \K 1A )am0 K + \K 1B }cos0 K , (124) 
|Ki(1400)) = \K 1A )cose K -\K 1B )sm9 K . (125) 

In flavor SU(3) symmetry limit, these mesons can not mix with each other; but since s quark is 
heavier than the u, d quarks, /^i(1270) and i^i(1400) are not purely l 3 Pi or l x Pi states. In general, 
the mixing angle can be determined by experimental data. The partial decay rate for r~ — ► K\v T 
is given by: 

r(T - _ KlKr) = £* rVrM (1 - §y (1 + m , (120) 
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with the measured results for branching fractions 

BR(t~ -> Xi(1270)i/ T ) = (4.7 ± 1.1) x 10" 3 , BK{t~ -> Ki(1400> r ) = (1.7 ± 2.6) x 1(T 3 .(127) 
We can straightforward obtain the longitudinal decay constant (in MeV): 

1/^(1270)1 = 169^?; 1/^(1400)1 = 1251^. (128) 

In principle, one can combine the decay constants for Kia, K\b evaluated in QCD sum rules with 
the above results to determine the mixing angle 9k- But since there are large uncertainties in 
Eq. (1128j) . the constraint on the mixing angle is expected to be rather smooth: 

- 143° < 6 K < -120°, or - 49° < K < -27°, or 37° < 9 K < 60°, or 131° < K < 153129) 

where we have taken the uncertainties from the branching ratios in Eq. (|127p and the first Gegen- 
bauer moment af" into account but neglected the mass differences as usual. For simplicity, we use 
two reference values in Ref. 

9 K = ±45°. (130) 

Besides, the flavor-octet and the flavor-singlet can also mix with each other: 

|/i(1285)> = |/i)cos03 Pl + |/ 8 }sin03 Pl , |/i(1420)) = -|/i)sin03 Pl + |/ 8 )cos03 Pl , (131) 
|/ii(1170)) = |/ii)cos6»i Pl + |/t 8 )sin0i Pl , |/ii(1380)) = -\h 1 )sm9i Pl + |/i 8 )cos0i Pl . (132) 
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The references points are chosen as: 63 Pl = 38° or 63 Pl = 50°; 6i Pl = 10° or 61 Pl = 45° [37j]. We 



should point out that if the mixing angle is 9 = 35.3°, the mixing is ideal: /i(1285) is made up of 
uu+dd wj^ig ^(1420) is composed of ss. Thus some of the form factors are very small which will 
of course give small production rates of this meson. 

Apart from these differences, the expression for B — > is different from B — > V"f in more 
aspects. Since the twist- 2 LCDA <j>\\ is normalized to ajj, we should replace the decay constant 
fv by Iao^q in the first two annihilation diagrams. As there is no overlap between an axial-vector 
meson and a photon (wrong parity), there is no contribution from the two photon electromagnetic 
operator diagrams in Fig. [8J Regardless of these differences, the factorization formulae for B — ► 
decays can be obtained from the corresponding B — > V7 ones using the replacement in Eq. (|39|) if 
the electroweak current is cr^ u (l + 75) or 7^(1 — 75) type. If the current is a^ u (l — 75), we should 
add an additional minus sign. In annihilation diagrams, if the electroweak current is LL or SP in 
the lower two diagrams, we need replace the distribution amplitudes as in Eq. ()39j) : while we add 
a minus sign if the current is LR. We show the formulas in Appendix ICl 

Branching ratios (in unit of 10 -6 ) and direct CP asymmetries (in %) for B — ► (a±, 61)7 processes 
are calculated straightforward as follows: 

(133) 
(134) 
(135) 
(136) 
(137) 
(138) 

i4&(B--ar(1260) 7 ) = II^^H^ls, (139) 
A% P (B° - a?(1260) 7 ) = S.Sl^t™;™, (140) 
Af P {B 8 - a?(1260) 7 ) = O.Si^i?;^;™ (141) 

A& P {B- - &r (1235)7) = 16.0iSttlKI^. (I 42 ) 
A d SHB° - 6?(1235) 7 ) = H-Oig^H;^, (143) 
A&(B, - 65(1236)7) = -0.5lg^£g£H, (144) 

while we give branching ratios and CP asymmetries for B — > Ki(f\, ^1)7 in table ITX1 and IX l respec- 
tively. The errors for these entries correspond to uncertainties in the input hadronic quantities, 



BK(B~ - 


> (1260)7) 


n +l. 6+0.4+0.4+0.8 
— o.U_ 1 _ 1 _ _3_Q >2 _o.7) 


BTZ(B° - 


-> o§ (1260)7) 


1 c+0.7+0. 2+0.2+0.4 
~~ i - d -0.5-0.2-0.1-0.4' 


BK(B S - 


-> o? (1260)7) 


_ ( i+0. 6+0.3+0. 1+O.On v 1f1 -4 
— l z -- L -0.5-0.1-0.1-0.oJ x 1U 


BK{B~ - 


*&r (1235)7) 


9 n +l. 0+0.4+0.3+0.6 
— 4.U_o.7_0.3-0.1-0.5> 


BTZ(B° - 


- 6?(1235) 7 ) 


1 1+0.5+0.2+0.2+0.3 
~~ J -- L -0.3-0.1-0.1-0.2' 


B1Z{B S - 


- 6?(1235) 7 ) 


— /K 4+1.0+6.4+0.3+2.1\ v -i n -5 
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TABLE IX: CP-averaged branching ratios (xlO -6 ) of B — > A7 decays obtained in pQCD approach using 
two different mixing angles; the errors for these entries correspond to uncertainties in the input hadronic 
quantities, from the scale-dependence, and CKM matrix elements, the Gegenbauer moments of the axial- 
vector mesons respectively. 



Modes 


K = 45° 


0k= -45° 


Exp. 


B~ -> iff (1270)7 


1 c 1 n 1 1 a 1 /ii 

1 "ZA +68+21+4+41 
io4: -49-18-4-38 


1 ,1 +i.z+U.o+U.0+5.U 
1 * ft — 0. 7-0. 6-0. 0-2.0 


42.8 ±9.4 ±4.3 [21] 


B° -> ^(1270)7 


1/11 +64+19+4+45 
1^1-48-18-4-41 


1 r\ O 1 O O 1 M O 1 E A 

-1 y| +u.y+u.o+u.u+o.4 

0.6-0.5— 0.0-1.9 




B~ -> iff (1400)7 


1 A -\- 1 . atU . OtU . UtO . u 

1 -*-0. 7-0. 6-0. 0-2.0 


1 Q/1+ D ' 5 + ' il T^T-il 

10 ^-49-18-4-38 


< 14.4 


B y) — > if ^(1400)7 


1 ,+0.9+0.3+0.0+5.4 
^-O. 6-0. 5-0. 0-1. 9 


I/it 4-644-19+4+45 
1^1-48-18-4-41 




B s AJ ) (1270)7 


... q+O. 07+0. 02+0. 02+0. 34 
U - ly -0. 06-0. 03-0. 01-0. 22 


- oo+0. 24+0. 09+0. 07+0. 44 
"•^"-0. 15-0.07-0.03-0.32 




B s -> if? (1400)7 


pi qq+O. 24+0. 09+0. 07+0. 44 
U - JS -0. 15-0.07-0.03-0.32 


„ . q+0. 07+0. 02+0. 02+0. 34 
U - ly -0. 06-0. 03-0. 01-0. 22 




Modes 


6»3 Pl = 38° 


03 Pl = 50° 




B° -> /x (1285)7 


-1 7+U.8+U.2 + U.2+U.0 

L - 1 -0.6-0.2-0.1-0.4 


1 fi+U. /+U.2+U.Z + 0.4 

x -°-0. 5-0.2-0. 1-0.4 




B° -> /1 (1420)7 


( A q+2.3+0.6+0.7+3.9n 1 n -3 
l 4 - y -l. 7-1. 1-0.3-2. 7J x 1U 


n 1 1 +0.05+0.01+0.02+0.04 
u - 11 -0.04-0.02-0.01-0.04 




5° -> /1 (1285)7 


n -,1+0.05+0.01+0.00+0.03 
u - 11 -0.04-0. 01-0. 00-0. 03 


q a+1.6+0.4+0.1+0.7 
°-°-l. 2-0. 4-0. 1-0.7 




A (1420)7 


R i q+24.5+5. 5+1. 8+17.4 
oi - a - 18.9-6.0-1.8-15.5 


ro 9+22.9+5.1 + 1.6+16.7 
JO - z - 17. 7-5. 6-1. 7-14. 8 




Modes 


6i Pl = 10° 


9i Pl = 45° 




B° -> fc x (1170)7 


n QQ+0.43+0. 16+0. 14+0.24 
u - yJ -0. 33-0. 13-0. 06-0. 21 


1 9 ,+0.55+0. 20+0. 18+0. 31 
-■-■ zq: -0. 41-0. 16-0. 08-0. 27 




5° -> ^(1380)7 


n 9O+0-12+0. 05+0. 04+0. 07 
u - zo -0. 09-0. 04-0. 02-0. 06 


/ n +0.8+0.3+0.3+0.3\ v 1n -2 
l z - u -0.7-0.3-0.1-0.3i x lu 




B° -> ft x (1170)7 


7 q+2. 9+1. 0+0. 2+1. 8 
' - y -2. 2-0. 7-0. 2-1. 6 


9 q+0. 9+0. 3+0. 1+0. 7 
z - J -0.7-0.3-0.1-0.6 




5° -> /ii (1380)7 


, + 16.8+5.6+1.3+11.0 
**'*-12. 8-4.1-1.3-9. 7 


m n + 18.8+6.3+1.5+12.2 
ou - u - 14.3-4.5-1.5-10.7 





the scale-dependence, CKM matrix elements, and the Gegenbauer moments of the axial-vector 
mesons. It is noted that theoretical uncertainties for branching ratios are quite large. The branch- 
ing fractions of B — ► ai(1260)(6i (1235))7 are larger than that of B — > pj, as we have shown that 
B — > A form factors are larger. As we have mentioned in the above, there are some ambiguities in 
the quark content of B — ► Ki(fi,h\)^: these mesons are mixtures but the mixing angles are not 
uniquely determined. The reference points for the mixing angles are two-fold, thus we give two 
different kinds of results collected in these two tables. The branching ratios and CP asymmetries 
are very sensitive to the mixing angles, which is not quite constrained. 



Experimentalist gave results for B~ — ► iff (1270)7 [211 ] shown also in table HXl Compared with 
it, our result for B~ — > iff (1270)7, i s about 3 times larger, when Ok = 45°; or very smaller than the 
experimental results, when Ok = —45°. In Fig. [91 we show the strong dependence of the branching 
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TABLE X: Direct CP asymmetries of B — > A-f decays obtained in pQCD approach using two different 
mixing angles; the errors for these entries correspond to uncertainties in the input hadronic quantities, from 
the scale-dependence, and CKM matrix elements, respectively. 



Modes 


0jf - 45° 


9 K = -45° 


B~ -> iff (1270)7 


-.inninoinnini 
fi + U.U-f-U.-i-r-U.U+U.l 

u ' u -0. 0-0. 1-0.0-0.1 


_ _ i n i -I Q_\_A n 

°-°-l. 3-0. 3-0. 3-13.0 


B° -> A^(1270) 7 


-0.2 ±0.0 ±0.0 ±0.0 ±0.0 


i n n i n i in n i n o 

u *- L -0.1-0.3-0.0-0.4 


B~ -> Jiff (1400)7 


i n "7_i_i ill o_\_A fl 
O Q + U. (+1.1 + 1.8-1-4.1) 

3-0.3-0. 3-13.0 


n n_i_n oj-f) n_i_n i 
"•"-0.0-0.1-0.0-0.1 


B° — > JK"J (1400)7 


_i_n 0-1-0 1 -|~n 0-1-2 9 
"^-0. 1-0. 3-0. 0-0. 4 


-0.2 ± 0.0 ± 0.0 ± 0.0 ± 0.0 


S s -> Af(1270)7 


o 4+0. 0+0. 5+0. 4+5.1 
°- 4 -2.9-3.4-0.4-11.2 


q 1+3.9+4.0+0.8+8.6 
• J ' 1 — 0. 1-3.1-0. 6-14.1 


S 5 -> K\\ 1400)7 


o 1+3.9+4.0+0.8+8.6 
^•-t-0. 1-3. 1-0.6-14.1 


a a +0.0+0.5+0.4+5.1 
~ a - 4 -2.9-3.4-0.4-11.2 


Modes 


03 Pl = 38° 


03 Pl = 50° 


B ^ A (1285)7 


i n a i n Q_i_n n i 
°-^-0. 1-0.5-0.2-0.1 


q A +U. ( -|-U.o-)-U. Z-\-V.o 

°'^-0. 1-0.4-0.2-0.1 


B° -► A (1420)7 


— - _j o ri-un T-i-fi ^-i-fi n 

'• 1 -3.1-7.7-0.5-2.3 


- L -0.3-1.9-0.3-0.4 


B° -> A (1285)7 


n i +0.3+0.1+0.0+0.1 
u - ± -0.0-0.1-0. 0-0.2 


n o+O.l+O.O+O.O+O.O 
u - z -0. 0-0. 1-0.0-0.0 


Bj-f/i (1420)7 


-0.2 ±0.0 ±0.0 ±0.0 ±0.0 


-0.2 ±0.0 ±0.0 ±0.0 ±0.0 


Modes 


0i Pl = 10° 


0i Pl = 45° 


5° -> /i! (1170)7 


in O+0.0+1.4+0.4+0.0 
iu - z -0.9-2.5-0.7-0.4 


in 1+0.1 + 1.7+0.4+0.2 
lu ' 1 -0.5-2. 3-0. 7-0. 3 


5° -> h x (1380)7 


q o+0. 8+2. 4+0. 4+1.1 
J -°-0. 0-2. 0-0. 7-0.0 


i i q+0. 0+0. 0+0. 5+0.0 
1J --°-5. 1-4.2-0.7-3.5 


B° s h!(1170)7 


-0.2 ±0.0 ±0.0 ±0.0 ±0.0 


n i +0.0+0.0+0.0+0.0 
u - J --0. 0-0. 1-0. 0-0.0 


B° -► /i! (1380)7 


-0.2 ±0.0 ±0.0 ±0.0 ±0.0 


-0.2 ±0.0 ±0.0 ±0.0 ±0.0 



ratio on the mixing angle. At Ok = 45°, the B~ — > K± (1270)7 receives almost a maximal 
branching ratio. The current B~ — > XT' (1270)7 experiment implies our chosen two reference 
points of the mixing angle are not favored. From this Fig. [9l we could read out the experimental 
constrained mixing angle value, which are also two fold. However large hadronic uncertainties and 



the missing next-to-leading order corrections 



801 ] plus the still large experimental error bars make 



this constraint not very effective. Except for these two processes, other decay modes have not 
been measured. Here we refrain from a direct comparison with the previous studies on B — > Aj 
?], Q, 81; 82, 83, 84], as the analysis is similar in B — ► A form factors which has been performed 



in section IIIII 
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FIG. 9: The 9k dependence of the B — > K 1 (1270)7 branching ratio. The region between the two horizontal 
lines are allowed by the experimental ler bound, where we add the statistic and the systematic uncertainties 
linearly: 29.1 < BTZ < 56.5. 

VI. SUMMARY 



pQCD approach is based on hr factorization where we keep the transverse momentum of valence 
quarks in the meson, to smear the endpoint singularity, fey resummation of double logarithms 
results in the Sudakov factor. Resummation of double logarithms from the threshold region leads 
to the jet function. Sudakov factor and jet function can suppress the contribution from the large b 
region and small x region, respectively. This makes the p QCD approach self-consistent. Inspired 
by the success of pQCD approach in non-leptonic B decays [85j], we give a comprehensive study on 
the charmless Br s \ — > 1^(^4)7 decays in pQCD approach. 

Semi-leptonic and radiative decays are somewhat simpler than non-leptonic decays as only 
one hadronic meson involved in the final state. In this case, the dominant amplitude can be 
parameterized in form factors. In order to make precise prediction and extract CKM matrix 
elements, we have to know the behavior of form factors. In pQCD approach, the final state 
meson moves nearly on the light-cone and a hard-gluon-exchange is required. Thus the dominant 
contribution is from the hard region which can be factorized. In section IIII) we have used the 
same input hadronic parameters with Ref. [n| and updated all the B — > V decay form factors in 
pQCD approach. Compared with the results evaluated from other approaches, we find: despite of 
a number of theoretical differences in different approaches, all the numerical results of the form 
factors are surprisingly consistent with each other. 

The 10 B — > V'j decay channels can be divided into three categories based on their dominant 
quark transition b — > s, b — > d and the annihilation topology. Our prediction on the first category 
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of decays BTZ(B — > if* 7) is consistent with the averaged value from experiments. On the b — > d 
transition, branching ratios have been given by BaBar and Belle collaborations with still large error 
bars. We find our results are well consistent with BaBar's results but a little far from the Belle's 
central value results in some channels. We also give our predictions on the purely annihilation type 
decays with very small branching ratios in SM. In three b — > s transitions B~ — ► K*~j, B° — ► K*°j 
and B s — > cpj, the direct CP asymmetry is small, since the tree contribution is suppressed by the 
CKM matrix element. In the b — ► d transitions B~ — > p~7, B° — > jf(uj)^/ and B s — > K* 0r y } 
the tree contribution can be comparable with the penguin contribution. Thus we obtain large CP 
asymmetries in these four processes. In SM, the two quantities Sf and Hf in time-dependent decay 
are expected to be rather small. This is due to the fact that the dominant contribution to decay 
amplitudes comes from the chiral-odd dipole operator O7. Except for a few decays discussed in the 
above, all other decay modes including B s decays and annihilation type decays, remain essentially 
unexplored. We wish a wealth of measurements at the B factories and other experiments in the 
future. 

In section we also study B — > A form factors. As the quark contents (to be precise the 
mixing angle) for the axial-vectors have not been uniquely determined, we give two different kinds 
of results for the form factors according to different mixing angles. For the axial-vector mesons 
/1, we have used the mixing angle between the octet and singlet: 9 = 38° (50°) which is close to 
the ideal mixing angle 8 = 35.3°. With this mixing angle, one can easily check that the lighter 
meson /i(1285) is made almost up of uu +^ d while the heavier meson /i(1420) is composed of ss. 
Thus partial decay widths of B — ► f\ (1420)7 an d B s — > f±( 1285)7 are suppressed by the flavor 
structure. In Fig. [U we show the strong dependence of the B~ — > iff (1270)7 decay branching 
ratio on the mixing angle 9k- Our calculation can be used to constrain this mixing angle using 
experimental measurements provided with well understood hadronic inputs. The study of higher 
resonance production in B decays can help us to uncover the mysterious structure of these excited 
states. 
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APPENDIX A: PQCD FUNCTIONS 

In this appendix, we group the functions which appear in the factorization formulae. The hard 
scales are chosen as 

t a = maxj-y/iJrnB, 1/bi, l/b 2 }, t' a = max{^/i7m B , l/&i, l/b 2 ], (Al) 



t b = max{^xix 2 m B , V i 1 ~ x 2)m B , 1/&1, 1/M> (A2) 



t' 6 = max{ x /xi2;2mB, a/1 + xim B , 1/&i, I/&2}, (A3) 
t c = max{ v / |xi - x 2 |m B , y/x2~m B , 1/61,1/62}, (A4) 
i' c = max{ v / |xi - x 2 \m B , ^fx{m B , l/&i, I/&2}, (A5) 



i d = max{^/i7x2"m B , I/61}, (A6) 

i e = max{Vl + x\m B , l/&i}, t' e = max{y / x7m£, l/&i}, (A7) 
i/ = max{Vl - x 2 m B , I/&2}, = max{y / i 2 "m_B, I/62}. (A8) 

The functions /ij in decay amplitudes are from the propagators of virtual quark and gluon and 
are defined by: 



d(B)K (VBm B b 1 )I (VBm B b 2 ) 
+e(-B)i^H^\V^Bm B b l )J (V^Bm B b 2 )\ + (61 ~ & 2 )}, (A9) 

where h£\z) = J„(z) + i Y (z). 

The Sudakov factor from threshold resummation is universal, independent of flavors of internal 
quarks, twists, and the specific processes. To simplify the analysis, the following parametrization 
has been used [261 ]: 



St(x) = 2l+2 Z^/l + c) [x(l - x)f , (A10) 

(1 + Cj 

with c = 0.4. This parametrization, symmetric under the interchange of x and 1 — x, is convenient 
for evaluation of the amplitudes. It is obvious that the threshold resummation modifies the end- 
point behavior of the meson distribution amplitudes, rendering them vanish faster at x — > 0. 
The evolution factors E^ and E$ are given by 

E e (t) = a s (t)S t (x 2 ) exp[S B (t) - S 2 (t)], E' e {t) = a s (t)S t ( Xl ) exp[-S B (t) - S 2 (t)\, (All) 
E a (t) = St(zi)exp[-S B (t)], E' a (t) = S t (x 2 )exp[-S 2 (t)}, (A12) 
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in which the Sudakov exponents are defined as 



m B \ 5 r dp 

' Xl ~^ J + 3 A b y 7 ^° s ^^' 



S 2 (t) = s[x 2 ^,b^+s^l-x 2 )^,b 2 ^+2j' b ^ q (a s (p)), 



(A13) 
(A14) 



with the quark anomalous dimension 7 g = —a s /7r. The explicit form for the function s(Q, b) is: 



s(Q,b) 



— ^-9 In 

2/?i \b 



AW 



q-b) + 



A® 



m \b 



ln(2g) + 1 ln(26) + 1 



+ 



8/3? 



A m A d) /e 2 7E -i 

— In 

4/3? 4/3i 

In 2 (2<7) - In 2 (26) 



where the variables are defined by 



In 



by 

(A15) 



q = ln[Q /(V2A)], 6 = In[l/(6A)], 



(A16) 



and the coefficients A® and /3j are 



33-2 ra/ 153 -19n, 

Pi = tt; — > P2 = 



12 

2 



24 



, m 4 , f2l 67 tH 10 8 _ . .1 
= -, A (2 > = n f + -/3iln(- 



3 27 



(A17) 



ra^ is the number of the quark flavors and 7_b is the Euler constant. We will use the one-loop 
running coupling constant, i.e. we pick up only the four terms in the first line of the expression for 
the function s{Q,b). 



APPENDIX B: ANALYTIC FORMULAE FOR THE B -► V7 DECAY AMPLITUDES 

The analytic formulae for B~ — ► p~7 is: 

A\B~ - P ~i) = ^Lv^V^M^ + Mt\Qu) + M\ u + M± n LL H ai ,Q u ) + A^^K Q d , Q u )} 
+ ^V cb V c * d {M\^ + Mf(Q u ) + Ml c } 

-^=V tb V t * d {M^ + Mlf + -Mg } (Q u ) + A^i a ^ L) («4 + aio, Qu) 
+Mlf>(«4 + aio, Qd, Qu) + M^r)(a 6 + a 8 , Q d , Q u )}, (Bl) 



44 



while the expression for B — > if* 7 is basically the same except with the only difference in the 
CKM matrix elements: V qd — > V qs . The formulas for other channels are 

V2A\S° - A) = ^V ub V: d {M± n LL \a 2 , Q d ) + Ad 6 n LL) («2, Qu, Qu) - - M^(Q d ) - M\ u ) 

+ ^V cb v: d { - Jrfg - Mf(Q d ) - M 2c ) - ^V tb V t * d { - Af 77 - Mlf 

-Mlf(Q d ) + M± k LL H~a, + \a 7 + ^a 9 + i Ol0 , Q d ) 

(a 3 + 09, Qu, Qu) + M±£ R) (a s + a 7 , Q u , Q u ) 
+M a ( ^ L) (-a 3 - 04 + ^a 9 + -010, Qd, Qd) + M a ^ R) (-a 5 + ^a 7 , Q d , Q d ) 
+-^l S n P) (-«6 + \a % , Q d , Q d ) + MUQu -Qd)}, (B2) 



V2A l (B° - uyy) = -^V ub V: d {M± k LL \a 2 , Q d ) + M± n LL \a 2 , Q u , Q u ) + + M\^(Q d ) + M 2u ) 

+^V cb v: d {M\ { : ) + M^(Q d ) + M\ c ) - ^V tb V* d {M 7l + 

+M i ff(Q d ) + M±f L) (2a 3 + a 4 + 2a 5 + X -a 7 + ^a 9 - -a 10 , Q d ) 

+M2'n L) (a 3 + a 9 , Q u , Q u ) + Mf n > n LR \a 5 + a 7 , Q u , Q u ) 

+-^l ( nn LL) («3 + a 4 - ^9 - \a lQ , Q d , Q d ) + ^i^ Li?) («5 - \a 7 , Qd, Qd) 

+MW?(ae - ^as, Q d , Qd) + MUQu + Qd)}, (B3) 



^(S°-r°7) = %V ub V: s {M i i: ) +M^(Q d )+M\ u } 



Gf 



^v cb v: s {M i i: ) +M i ic ) (Qd)+M 2c ) 



-^L Vtb V t * s {M 7j + + Mlf(Q d ) + M±£ L \a A - \a w , Q d ) 

+M'tn L) (a A - ±a w ,Q s ,Qd)+M2P(a 6 - \a 8 ,Q s ,Q d )}. (B4) 

The expression for B s — > if *°7 can be obtained by replacing V qs by from 5° — > K*°j. 
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The formulas for the B s — ► ^7 decay are 



+^v cb v: s {u i i: ) +Mt\Q s )+M^} 

-^V tb VZ{M\ 7 +Mif+Mf(Q s ) 

+^nn L) (a 3 + a A + a 5 - X -a 7 - ^a 9 - -a w , Q s ) 

+M i ^ L) (a 3 + a 4 - \a 9 - -a 10 , Q s , Q s ) + M±£ R) (a 5 - -a 7 , Q s , Q s ) 

+M^\ae - \a 8 ,Qs,Qs) +M en {Q s )}. (B5) 

For the annihilation type decays, we have 

A\B ^ch) = -^V tb V; d {M^ LL \a 3 + a 5 -^a 7 -^a 9 ,Q s )+Ml n (Q s ) 

+Mtn LL \as ~ la 9 ,Q s ,Q s )+Mtn LR) (a 5 ~ \a 7 ,Q s ,Q s )}, (B6) 



V2A\B s ^p^) = ^V ub V: s [M^ LL \a 2 ,Q s )+M^ L \a 2 ,Q u ,Q u )) 



V2 

^V tb V*{Mt k LL H^a 7 + °-a 9 , Q s ) + M en {Q u - Q 

+-^l b n LL) («3 + 09, Qu, Qu) + Mt'n R H^ + «7, Qu, Qu) 

+M±' n LL H-a 3 + ±a 9 ,Q d ,Q d )+M±' n LR \-a 5 + ^a 7 ,Q d ,Q d )}, (B7) 

V2A i (B s — ► ury) = ^V ub V: s {M± k LL Ha 2 ,Q s ) + Mlf n LL \a 2 ,Q u ,Q u )} 

-^V tb V t * s {M± k LL \2a 3 + 2a 5 + ^a 7 + ^ag, Q s ) + Af en (Q u + Q d ) 

+Mtn LL) (a 3 + a 9 , Q n , Q u ) + M i ^ R \a b + a 7 , Q u , Q u ) 

+M l ± LL \a 3 - \a 9 ,Q d ,Q d )+M^ n LR \a 5 - \a 7 ,Q d ,Q d )). (B8) 

APPENDIX C: ANALYTIC FORMULAE FOR THE B -> A7 DECAY AMPLITUDES 

The expression for B — > l Pi 7 is different from P — > V7 in various aspects. The first two 
annihilation diagrams vanish for neutral axial vector mesons (not including K\b)- There is not any 
two-photon diagram contribution in B — > decays. The explicit formula for the P~ — > 6^ (1235)7 
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decay amplitude is 



A\B- - 6^(1235)7) = ^V ub V: d {M\ { : ] +Mt\Q u )+Ml u + M'± LL \ ai ,Q d ,Q u )} 

+^VM{Mf? + Mf(Q u ) + Ml c } 
-^V tb V t * d {M^ + Mlf + ^««(0 U ) 
+Mt U) (a4 + «io, Q d , Qu) + M l ^P(a 6 + a s , Q d , Q B )}, (CI) 

while the expression for U-spin related process B~ — > K^ B ^ is basically the same except with the 
only difference in the CKMmatrix elements: V qd — > V qs . The formulas for the neutral decays modes 
are 

V2A\B° - 6°(1235) 7 ) = ^VM{m$£ L ' '(a 2 , Q u , Q u ) - - MfJ{Q d ) - M\ u ) 



V2 

-M^\Q d ) + Mtn L) ^3 + a 9 , Q u , Q u ) + M± n LR \a 5 + a 7 , Q u , Q u ) 
+-^l ( nn LL) (-a 3 - a 4 + \a 9 + -a 10 , Qd, Qd) - M 



i{a) 



+M±£ R H-a 5 + \a 7 ,Q d ,Q d )+M a ^P(-a 6 + \a % ,Q d ,Q d )}, (C2) 



A\B° - K&y) = ^Z Vub v: s {M[M + M^(Q d ) + M\ u ) 
+^V cb v: s {Mt ) +Mf\Qd)+M\ c ) 

-^=V tb V t *{M 77 + Mlf + M^(Q d ) + Mti: L \a, - ±a 10) Qd) 
+M%£ l H<h - \a w ,Q s ,Q d ) + M a ^{a % - ±a 8 ,Q s ,Q d )}, (C3) 

and the expression for B s — ► K® B j can be obtained by replacing V qs by V qd . 
The decay amplitudes involving /ii(1170) and /ii(1380) can obtained as 

A\B -> /ii (1170)7) = A*(.B -> /ii7)cos0i Pl + A\B -► /i 8 7)sin0i Pl , (C4) 
A* (5 -> /ii (1380)7) = -A*(B fri7)sin0i Pl + A(B fr 8 7)cos0i Pl , (C5) 
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where B denotes B° or B® with 

V2A\B°^h 8l ) = ^VM^Mtn^Ha^Q^Q^+M^+MtHQ^+M^} 

-^V tb V t * d -L{M^ + M*ff + Mi s ;'(a 6 - la 8 , Qd, Qd)} 
+M i $(Q d ) + M i ^ L \a 3 + afl, Q u , Q u ) + M^ R) (a 5 + a 7 , Q U! Q u ) 
+-^l ( nn L) («3 + a 4 - la 9 - -a 10 , Q d , Q d ) + M a ^ R) (a s - la 7 ,Q d ,Q d ) 
-2MfJ L \a 3 - l -a 9 ,Q s ,Q s ) - 2M a ^ LR \a 5 - l -a 7 ,Q s ,Q s )}, (C6) 

y/2A*(2°^hn) = ^V ub v: d ^{Mtn LL Ha 2 ,Q u ,Q u )+M\ i : ) +M^(Q d )+M^ 
+^VM±{M l l:~> + Mf(Qd)+M 2c } 



G F 



V tb V t * d ^={M 7l + + M a ^P( a(i - ±a s , Q d , Q d )} 



V2 10 td Vs 

+M i ff(Q d ) + M a ^ n LL \a 3 + a 9 , Q u , Q u ) + Mf n > n LR \a 5 + a 7 , Q u , Q u ) 
+-^l ( nn LL) (a 3 + a 4 - ^a 9 - -aio, Qd, Qd) + Mf L ^ R) {a 5 - -a 7 , Q d , Q d ) 
+Mt'n L \a 3 - lag,Q s ,Q s )+Mt n LR \a 5 - \a 7 ,Q s ,Q s )}, (C7) 



A\B s ^h %1 ) = ^L Vub V: s ^{M\ { : ) +M\^(Q s )+M 2u } 



^V cb V: s MM i i: ) +M i ic ) (Qs)+M 2c } 



V2 u ° us Vs 

-^V»VZ=£{M^ + Mlf + M^(Q S ) + MfrVe - \as, Q s , Q s ) 
+Mf n ' n LL \a 3 + a 4 - lag - la w , Q s , Q s ) + M^ R \a^ - \a 7 ,Q s ,Q s )], 



+ ^V ub V: s ^{Mtn L) (a 2 ,Q u ,Q u )} 

-^V tb V*^{M a f n LL \a 3 + ag, Q u , Q u ) + MfJ R \a 5 + a 7 , Q u , Q u ) 
+Mt'n L) (a 3 - lag,Q d ,Qd)+M±' n LR \a 5 - \a 7 ,Q d ,Qd)}, (C8) 
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AiB.^hn) = ^v ub v: s ^={M\ { : ] +M\ { *\q s )+mQ 

+^V cb V c * s ^={M^ + Mf(Q s )+Ml c } 

-9z Vtb V*-L{M^ + Mlf + Mlf(Q s ) + M2P(a, - \a 8 , Q„ Q s ) 
+Mt' n LL \a 3 + a, - l -a 9 - \a 1Q ,Q s ,Q s ) + M%£ R Ha 5 - \a 7 ,Q s ,Q s )}, 
+ ^V ub V: s ^={Mtn LL) (a2, Qu, Qu)} 

-^V tb V t * s ^={Mt'n L) (a 3 + 09, Q u , Q u ) + M^ R) (a 5 + a 7 , Q u , Q u ) 
-Mf n ' n LL \a 3 - ±a 9 ,Q d ,Q d )+M±' n LR \a 5 - \a 7 ,Q d ,Q d )}. (C9) 



The annihilation type decay amplitude is: 



y/2A\B a - &?(1235) 7 ) = ^KfcKs^l^K Q u , Q u )} 



-^V tb V t * s {M±> n LL \a 3 + 09, Q u , Q„) + ^l 6 ^(a 5 + a 7 , Q u , Q tt ) 
+M l ^ L \-a 3 + ±a 9 , Q d , Q d ) + M^ R \-a^ + ±a r , Q d) Q d )](C10) 

In the following, we will give the analytic factorization formulae for i? — > 3 Pi 7 which is similar 
with B — > T^7 except for some differences in the flavor structure and zero contribution from the 
two-photon diagrams, 

^(p-^ar(1235) 7 ) = ^V ub V: d {Mt ) +Mt\Qu) + ML + M± k LL \a 1 ,Q u ) 



+Mtn LL \a 1: Q d , Qu)} + -^Vc^M^ + M\f(Q u ) + M\ c } 
-^VM{M\ 7 + Mlf + M^(Q U ) + M^ k LL \a 4 + a w , Q u ) 
+Mtn LL) (aA + ai 0j Qd, Qu) + -^a ( nn P) («6 + a 8 , Q d) Q u )}, (Cll) 
while the expression for P~ — ► -KT47 is the same except with the only difference in the CKM 
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matrix elements: V qd — > V qs . The formulas for other channels are 
sf2A\B° - a?(1235) 7 ) = ^V ub V: d {M^ LL \a 2 , Q d ) + M± n LL Ha 2 , Q u , Q u ) - 

-MT(Q d ) - M\ u ) + ^VM{ - - M*g\Q d ) - M 2c ) 

-M l 8 f(Q d ) + M± k LL \- ai + L 7 + L 9 + iaio, Q d ) 
+M a ^ L Ha 3 + 09, Q u , Qu) + ^± LR) (a5 + or, Qu, Qu) 
+-^l ( nn LL) (-a 3 - a 4 + ^a 9 + ia 10 , Qd, Qd) (C12) 
+-Mi 6 n LR) (-as + ^7, Qd, Qd) + A^i S „ P) (-a6 + ^a 8 , Q d , Qd)}, 



A\B^Kl Al ) = ^V ub v: s {M\^ + M^(Q d )+M 2u } 



V2 

+ ^V cb v: s {M\[ a) + Mf(Q d )+M 2c } 

-^ Vtb V*{M^ + Ml { ; ) +Mlf(Q d )+M± k LL \a, - \a w ,Q d ) 

+Mt£ L) (a 4 - ^a 10 ,Q s ,Q d )+M^\a 6 - ^a 8 ,Q s ,Q d )}, (C13) 

and the expression for B s — ► K^ A j can be obtained by replacing V qs by V qd . 
The decay amplitudes involving /i(1285) and /i(1420) can obtained as 

A\B -> /i (1285)7) = A7)cos03 Pl + / 87 )sin^ Pl , (C14) 

-> hi (1420)7) = -^(-B -► /i7)sin0 3Pl + A*(J3 -» f 8 j)cos9 3Pl , (C15) 

where -B denotes B° or with 

^(^-/st) = ^V^^M^Ha^Q^+Mt^^Q^Q^+M^ 

+Mt ) (Q d )+M l 2u } + ^LvM^M^+Mfm + Ml} 

+M±i l LL) (2a 3 + a 4 + 2a 5 + ^a 7 + ^a 9 - -oi , Qd) 

+-A< ( nn LL) («3 + a 9 , Qu, Qu) + (o 5 + a 7 , Qu, Qu) 

+-^l b n L) («3 + 04 " ^9 " 2«10, Qd, Qd) + -Mi ( l L/?) («5 " ^7, Qd, Qd) 

+M i ^ ) (a 6 - -a s , Q d , Qd) ~ 2M±f L) (a 3 + o 5 - ^a 7 - ^a 9 , Q s ) 
-2^i^ L )(a 3 - \a 9 ,Q s ,Q s ) - 2M i ^ R \a b - ^a 7 ,Q s ,Q s )}, (C16) 
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x/2^(S°^/ l7 ) = ^V ub V: d ^{M± k LL \a 2 ,Q d )+Mt n LL \a 2 ,Q Ul Q u )+M\ { : ] 

+M^(Q d )+Mi u } + ^VM-L{M[M + M l ^(Q d ) +M\ C ] 

+Mlf(Q d ) + M a ^ LL \2a 3 + a 4 + 2a 5 + ^a 7 + ^a 9 - ^a 10 , Q d ) 

+-^l ( nn LL) («3 + Ofl, Qu, Qu) + M a ^ LR \a 5 + a 7 , Q U) Q u ) 

+*K ( nn L) (a 3 + a 4 - ^ag - ^10, Qd, Qd) + -Mi ( l LR) («5 - \a 7 , Q d , Q d ) 
+M i ^P{a & - -a 8 , Q d , Qd) + M±^ LL) (a 3 + a 5 - ^a 7 - ^a 9 , Q s ) 
+Mt'n L \a 3 - \a 9 ,Q s ,Q s )+M±' n LR \a 5 - \a 7 ,Q s ,Q s )}, (C17) 

V2A l (B s ^f 8l ) = 9z Vub V: s -^={M± k LL \a 2 ,Q s ) + M a % n LL \a 2 ,Q u ,Q u )} 

-^V tb V*±{M a ^ LL H2a 3 + 2a 5 + \a 7 + \a 9 , Q s ) 
+Mt£ L) (a 3 + a 9 , Q u , Q u ) + M^ iR) («5 + a 7 , Q u , Qu) 
+ Mt'n L) (a 3 - 7j a 9 ,Q d ,Qd)+Mtn LR) (^ ~ \a 7 ,Q d ,Qd)} 
+^V ub v: s ^{Mt ) +M^(Q s ) + Mi u } 

+^V cb v: s ^{M\^ + M^(Q s )+Ml c } 

-^ Vtb V; s ^{M\ 1 +Ml { ; ) +M^\Q s ) + M2n\a 6 - \a,,Q s ,Q s ) 

+-^l a A LL) («3 + a 4 + a 5 - \a 7 - \ 9 - -a w , Q s ) (C18) 
+M±£ L) (a 3 + a 4 - ^a 9 - ^a 10 , Q s , Q s ) + M' a ^ n LR \a 5 - -a 7 , Q S ,Q S )}, 
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-^V tb V t *^={M±f L \2a 3 + 2a 5 + \a 7 + ±a 9 , Q s ) 
+M i ^ L \a 3 + a 9 , Q u , Q u ) + M^ R \a B + a 7 , Q u , Q u ) 
+Mt'n L) (a 3 - \a 9 ,Q d ,Q d )+M±' n LR \a 5 - \a 7l Q d ,Q d )} 

+^V ub V: s ^={Mt ) +M^(Q s )+Ml u } 

-^V tb v; s ^{M^+Mlf+Mlf(Q s ) + M^ n P) ^ ~ \^QsM 
+M±£ L \a 3 + a A + a 5 - ^a 7 - ^a 9 - -a w , Q s ) (C19) 
+Mt'n L) (a 3 + a 4 - ^a 9 - ±a 10 , Q s , Q s ) + -M^K - \a 7 , Q S ,Q S )}. 

For annihilation type decays, we have 

V2A i (B a — ► a? (1235)7) = ^K b C{^i a ^ L) («2,Qs) + ^i b n ii) (a 2 ,Q n ,Q u )} 

-^^{^^(^ar + ^a 9 ,Q s ) 

+Mi ( i Li) (a 3 + 09, Q„, Q«) + ^1^ LR) (a 5 + a 7 , Q u , Q tt ) (C20) 
+M* Li) (-«3 + ^a 9 ,Q d ,Q d ) +M^)(-a 5 + ±a 7 , Q d , Q d )}. 

Decay amplitudes of A*(B -> #1(1270)7) and #1(1400)7) ( nere 5 denotes # M)S and #1 

denotes (Kf)) can be obtained by: 

A\B -> #i(1270) 7 ) = sin(^)^(S #1^7) + 008(^)^(5 #i B7 ), (C21) 
A* (5 -> #1 (1400)7) = -*fa{0 K )A i (B^K 1 Bi)+cKx(0 K )A\B^KiA'i). (C22) 
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